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f~| . Abstract 

> . 

S In this paper we produce families of complete non compact Riemannian metrics with positive constant 
afc-curvature equal to 2 _fc ('^) by performing the connected sum of a finite number of given n-dimensional 
Delaunay type solutions, provided 2 < 2k < n. The problem is equivalent to solve a second order fully 
t-H ■ nonlinear elliptic equation. 
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O ■ 1 Introduction and statement of the result 

o : 

t-H . In recent years much attention has been given to the study of the Yamabc problem for a^. -curvature, 

briefly the <7/c -Yamabe problem. To introduce the analytical formulation, we first recall some background 

• *h . material from Ricmmanian geometry. Given (M, g), a compact Riemannian manifold of dimension n > 3, 

we denote respectively by Ric g , R g the Ricci tensor and the scalar curvature of (M,g). The Schouten 
^ . tensor of (M, g) is defined as follows 

A 9 ~ 1T=2 ( Rtc 9 - 2^=n R a a) ■ 

If we denote by Ai, . . . , X n the eigenvalues of the symmetric endomorphism g~ 1 A g , then the Ufe-curvature 
of (M,g) is defined as the fc-th symmetric elementary function of Ai, . . . , A„, namely 

o-k{g~ 1 Ag) := 2_, ' ••• ' for 1 < fc < n and ao(g~ 1 A g ) := 1. 

ii < ... < i fc 

The CTfe-Yamabe problem on (M, g) consists in finding metrics with constant -curvature in the same 
conformal class of g. The case fc = 1 is the well known Yamabe problem, whose progressive resolution is 
due to Yamabe [23], Trudinger [22], Aubin [1] and Schoen [19]. In order to present the existence results 
for fc > 2, when the equation becomes fully nonlinear, we recall the notion of fc-admissibility, which is a 
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sufficient condition to insure the cllipticity of the equation. A metric g on M is said to be fc-admissible 
if it belongs to the fc-th positive cone I^" , where 

geT+ a ] (g- 1 A g )>0 for j = l,...,k. 

Under the assumption that g is fc-admissible the (positive) Cfc-Yamabe problem on closed manifolds has 
been solved in the case k = 2, n — 4 by Chang, Gursky and Yang [5] [6], for locally conformally flat 
manifolds by Li and Li [13] (see also Guan and Wang [8]), and for 2k > n by Gursky and Viaclovsky 
[9]. For 2 < 2k < n the problem has been solved by Sheng, Trudinger and Wang [21] under the extra- 
hypothesis that the operator is variational. We point out that for k = 1,2 this hypothesis is always 
fulfilled, whereas for k > 3 it has been shown in [2] that this extra assumption is equivalent to the locally 
conformally flatness. 

Of interest in this paper is the construction of complete non compact locally conformally flat metrics 
with constant (positive) <7fc-curvaturc. These solutions can be regarded as singular solutions for the 
(Tfe-equation on the complement of a discrete set A on the standard n-dimcnsional sphere. To put our 
result in perspective, we recall that for k = 1, the first examples of conformal constant (positive) scalar 
curvature metrics with isolated singularities have been obtained by Schoen in [20]. Later, Mazzeo and 
Pacard (see [14] and [15]) have produced different families of solutions on the complement of a singular 
set A consisting of a finite disjoint union of closed smooth submanifolds of arbitrary dimension between 
and (n — 2)/2. Another existence result, in the case where A is given by an even number of points, is 
due to Mazzeo, Pollack and Uhlcnbcck [16]. We will return on this later, since our construction is closely 
related to their work. For 2 < k < n/2, the first examples of complete non compact metrics lying in the 
fc-th positive cone and having constant avcurvature have been obtained by the first author in a joint 
work with Ndiayc [17], assuming that the points of the singular set have a symmetric disposition. 

Some comments about the asymptotic behavior of the singular solutions are now in order. For k = 1 it 
follows from the works of Caffarelli, Gidas and Spruck [3] and Korevaar, Mazzeo, Pacard and Schoen [12] 
that every complete non compact locally conformally flat metric with constant positive scalar curvature 
must be asymptotic to a radial solution. In a recent work, Han, Li and Teixeira [10] have shown that this 
fact is true also for metrics of constant curvature lying in the fc-th positive cone, provided 2 < k < n/2. 
Notice that for k > n/2 the singularity is always removable. For these reason, it is clear that complete 
radial solutions arc going to play a fundamental role in our construction. These particular solutions are 
also known as Dclaunay-type metrics and have been classified by Chang, Han and Yang in [7] and we 
recall them briefly in Section 3. Essentially, they are conformally cylindrical metrics with a periodic 
conformal factor, whose minimum will be referred as Delaunay-parameter. 

As anticipated, the construction presented in this paper is inspired by [16] and consists in performing 
the connected sum of a finite number of Dclaunay-type metrics. The solutions obtained in this way are 
quite different from the ones produced in [17], which roughly speaking looks like a spherical central body 
with several Delaunay-type ends having small Delaunay parameters. In the present construction, the 
Delaunay parameters are not forced to be small, hence our solutions may possibly belong to a different 
connected component of the moduli space. 

To fix the notations, we recall that the connected sum of two n-dimensional Ricmannian manifolds 
(Di,gi) and (Z?2,<?2) is the topological operation which consists in removing an open ball from both 
Di and D 2 and identifying the leftover boundaries, obtaining a new manifold with possibly different 
topology. Formally, if pi <G Di and for a small enough e > we excise the ball B(pi, e) from Di, i = 1,2, 
the (pointwise) connected sum M e of D\ and Di along p± and pi with necksize e is the topological 
manifold defined as 

M e := D4 £ D 2 = [D 1 \%,e)UD 2 \% E )]/~, 

where ~ denotes the identification of the two boundaries dB(pi, e), i = 1, 2. Of course the new manifold 
M e can be endowed with both a differentiable structure and a metric structure, as it will be explicitly 
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done in Section 4. Even though from a topological point of view the value of the necksize is forgettable, 
it will be important to keep track of it, when we will deal with the metric structure. 

Concerning the solvability of the Yamabe equation (fc = 1) on the pointwise connected sum of manifolds 
with constant scalar curvature, we recall the results of Joyce [11] for the compact case and the already 
mentioned work of Mazzeo, Pollack and Uhlenbeck [16] for the non compact case. For 2 < k < n/2 and 
compact manifolds, a connected sum result has been provided by the first author in a joint work with 
Catino [4]. Our main result here is the following 

Theorem 1. Let (D\, gi), . . . , (-Dat, <7at) be a collection of n- dimensional Delaunay-type solutions (see 
Proposition 3.1) to the positive a^-Yamabe problem, with 2 < 2k < n. Then, there exists a positive 
real number e$ > only depending on n, k, and the C 2 -norm of the coefficients of gi and 52 such 
that, for every e € (0,£o], the connected sum M £ = D\$ £ . . .^ e Djq can be endowed with a metric ~g e with 
constant o^-curvature equal to 2~ fc ("). Moreover \\g e — gi\\c r (Ki) for any r > and any compact 
set Ki C Di \ {pi}, the Pi's, i = 1, . . . , N , being the points about which the connect sum is performed. 

Some comments about the strategy of the proof are in order. Incidentally, we notice that the constant 
^ _fc (fc) ar i ses naturally as the o^-curvature of the ro-dimensional standard sphere, so we will end up with 
a family of metrics {g e }e parametrized in terms of the necksize which satisfy 



To show the existence of these solutions, we start by writing down (see Section 4) an explicit family of 
approximate solution metrics {g £ }e (still parametrized by the necksize) on M £ . This metrics are complete 
and non compact, since they coincide with the original Delaunay-type metrics gi on Di \ B(pi, s), i = 1, 2, 
and are close to a model metric on the remaining piece of the connected sum manifold, which in the 
following will be referred as neck region. The metric which we are going to use as a model in the neck 
region is described in Section 3. It is a complete metric on M. x S"" 1 with zero cr^-curvature and yields 
a natural generalization of the scalar flat Schwarzschild metric. It has been successfully employed in [4] 
to treat the connected sum of constant scalar curvature manifolds and for the local analysis on the neck 
region we will refer to this work. 

The next step in our strategy amounts to look for a suitable correction of the approximate solutions to 
the desired exact solutions. This will be done by means of a global conformal perturbation. At the end it 
will turn out that for sufficiently small values of the parameter e such a correction can actually be found 
together with a very precise control on its size and this will ensure the smooth convergence of the new 
solutions g e to the former metrics gi on the compact subsets of Mi \ {pi}, i = 1, 2. We point out that 
it is also important to control the asymptotic behavior of such a perturbation in order to preserve the 
completeness of the approximate solutions. Typically, one is led to search for corrections which present 
a decay at infinity. 

The main point in the correction procedure is to provide invertibility for the linearized operator about 
the approxiamte solutions, together with uniform (with respect to the necksize parameter e) a priori 
bounds. This will enable us to carry out the perturbative nonlinear analysis (Section 7) by proving the 
convergence of a Newton iteration scheme. The uniformity of the a priori bound will follow from the 
use of weighted function spaces with a weighting function acting on the neck region, in analogy with the 
analysis contained in [4]. On the other hand, the invertibility issue is quite different from the compact 
case. In fact, in order to obtain the desired Fredholm properties for the linearized operator, we will further 
introduce weighting functions with gradient supported outside of a compact region of M e . In analogy with 
the case k = 1 (see [16]), the analysis is complicated by the lack of coercivity of the linearized operator. 
This is due to the conformal invariancc of the -equation. In fact, the functions which are responsible 
for this lack of coercivity arise as infinitesimal generators (Jacobi fields) of conformal transformations. 
As it will be made clear in Section 7, the geometrical interpretation of the Jacobi fields will be exploited 
in order to insure the completeness, after the perturbation, of the exact solutions. 
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2 Notations and preliminaries 

We fix now the notations that will be used throughout this paper. Let (M, g) be a compact smooth 
n-dimensional Riemannian manifold without boundary an let 2 < 2k < n. Taking advantage of this 
second assumption, we introduce the following formalism for the conformal change 

9u ■= g, 

where the conformal factor u > is a positive smooth function. In this context g will be referred as the 
background metric. At a first time the o~k -equation for the conformal factor u can be formulated as 

<9u l A- Bu ) = 2- k (l). 
We recall that the Schouten tensor of g u is related to the one of Ag by the conformal transformation law 

A 9u = A 3-^k u ~ 1 V 2 u+ jf2*jpu- 2 du®du- jJ^u-^dufg, 

where V 2 and | • | are computed with respect to the background metric g. For technical reasons, it is 
convenient to set 

D 9u ■ 2k u ^5" 

and to reformulate the afe-equation as 

(2-1) M{u,g) := a k (B- g J- (™)(i^) fe w ^ = 0. 

We notice that if two metrics g and g are related by g = (v j 'u) 4k : /( n ~ 2k ) g ; then the nonlinear operator 
enjoys the following conformal equivariance property 

(2.2) N(u,g) = (v/u)-£&M(v,g). 

The linearized operator of J\f( ■ , g) about u is defined as 

d 



(2.3) Uu,g)[w] :■- 

ds 



M (u + sw,g). 

s=0 



Most part of the analysis in this paper (Sections 5 and 6) is concerned with the study of the mapping 
properties of the linearized operator about the approximate solutions g e 's, that we will write in the form 
^4fc/(n 2k )g_ ^ s a dj rcc t consequence of the property (2.2), we have the following conformal equivariance 
property for the linearized operator 

(2.4) Hu,g)[w] = (v/u)-^Hv,g)[(v/u)w}. 
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3 Delaunay and Schwarzschild type metrics on R x S n 1 

We start this section with the description of a particular family of complete metrics on the cylinder 
R x §™ _1 with constant er^-curvature equal to 2 _fc ("). These metrics are conformal to the cylindrical 
one g cy i on the whole cylinder K x §™ _1 (notice that in the following, the cylindrical metric will also be 
denoted by g cy i = dt 2 + d6 2 , where dO 2 represents the standard metric on S n_1 ). 

Let us consider then on the standard cylinder (R x § n_1 ,di 2 + dO 2 ) a conformal metric g of the form 
g = w 4fe /(™-2 fc ) g cyli where the conformal factor v only depends on the t variable, i.e., v = v(t), and let 
us impose the condition at (g~ 1 A g ) = 2~ k (?) or equivalently 

(3.1) M(v,g cyl )=0. 

It is easy to observe that, under the usual change of coordinates, t = —\og\x\ and 8 = a;/|cc|, this 
corresponds to look for a metric on R™ \ {0} which has constant positive ovcurvature and which is 
radially symmetric. These metrics has been studied in [7] by Chang, Han and Yang and we refer the 
reader to their work for further details. Here we just recall the following 

Proposition 3.1 (Delaunay- type metrics). Let g v be a a metric onl" of the form g v = v 4k /( n ~ 2k ) g cy[j 
where v is a smooth positive function only depending on the variable t £ R. Let us define the quantity 

H(v,v) := [v 2 - (^)V] fc - v** 

Then, if H(v,v) = Hq E (0 , n 2 ^ 2 k ( "~ 2fc )" //2fc ), in correspondence of each Hq, there exists a unique 
solution v to 

(3-2) [^-{^kf^T-'iv-^fv] = ^t,**" 1 . 

satisfying the conditions v(0) = 0, and u(0) > 0. This family of solutions gives rise to a family of complete 
and periodic metrics on R x S"™ 1 satisfying 

(3.3) a k (B gv ) = 2- k Q inlxT 1 

This solution is periodic and it is such that < v (t) < 1 for all t £ R. In the following we will index 
the conformal factors and the metrics in this family by means of the parameter rj := v(Q) 2k ^ n ~ 2kS> which 

represents the neck-size. Notice that < r\ < ( "~ 2fc ) 1 ^ 2fc and that the period of vr> ri will be denoted by 
T v . 

As anticipated in the introduction, the first step in our strategy amounts to build approximate solutions 
on the connected sum of a finite number of Delaunay type solutions with possibly different neck-size 
parameter {D Vl ,g Vl ) 1 . . . , {D VN , g VN )- To this end, we need to modify the original metrics in a neighbor- 
hood of the points that we are going to excise, obtaining a new metric in the so called neck region. In 
the scalar curvature case (fe = 1), a clever choice turns out to be the (space-like) Schwarzschild metric. 
This is a complete scalar flat metric conformal to the cylindrical metric g cy i on R x § rl_1 . The explicit 
formula is given by 

4 

g := cosh(^t)"- 2 g cyl . 

In a similar way, it is easy to construct a complete conformal metric onRxS" -1 with zero o~k -curvature, 
for all 2 < 2k < n. We have 

Proposition 3.2 (Schwarzchild metrics). Let g v be a metric onRxS™ -1 of the form g v = v 4k ^ n ~ 2k ^ g cy i, 
where v is a positive smooth function depending only on t £E R. Let us define the quantity 

m ^v 2 (t)-(^.) 2 v 2 (t). 

Then, if h := h(0) > 0, the family of positive solutions to the equation 

a k (B gv ) = mlxT 1 

is given by v(t) = v^o cosh ( "~ 2fc t — c) , c G R. In the following, the solution with c = will be denoted 
by Vt.- 
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For sake of completeness and for future convenience, we recover from [4] the following formula for the 
linearized er^-operator about the Schwarzschild type metric. 



Lemma 3.3. The linearized Gk-operator about the o~k -Schwarzschild metric 

^°(vx,gcyl)[w] := — CTfc (B gg ) , 
as s=0 

where g s = g v +sw> is given by 

(3.4) h°(v^g cyl )[w] =-C n , k v s h k 1 T 1 [d? + 1 ^Ao- {^f" 



w . 



where C n , k = (^P)*. 

Incidentally, we note that the computation (see [4]) leading to (3.4) also shows that 

(3-5) v.-i-ABo) = {^) k - X - j 4"^ ¥ (*_?_,) ■ 

From this it follows that the (jfc-Schwarzschild metric <?£ belongs to H ri"_ 15 for 2 < 2k < n. 



4 Approximate solutions 

In this section we first describe the construction of the connected sum of a finite number of Delaunay-type 
solutions D ril , . . . , and then wc define on this new manifold a family of metrics which will represent 
the approximate solutions to our problem. 

Since the whole construction is local, we restrict ourself to the connected sum of two Delaunay type 
solutions (D^gi) and (D r)2 ,<72)- In the following we will denote by M e := D Vl %D V2 the manifold 
obtained by excising two geodesic balls of radius e € (0, 1) centered at p\ £ D m and P2 € D V2 and 
identifying the two left over boundaries. The manifold D m and D V2 are endowed with the metrics 

4fe 4k 

■9i = v D~™( dr i and 92 = v£~™ (drl + flsn-i) , 

respectively. Starting from g\ and 32, we will define on M £ a new metric g e which agrees with the old 
ones outside the balls of radius one around p\ and pi and which is modeled on (a scaled version of) the 
CTfc-Schwarzschild metric in the neck region. 

To describe the construction, we consider the diffeomorphisms given by the exponential maps 

exp pi : B(0 Pi ,l) c T Pi D m — > B{ Pll l) c D, n , i = 1,2. 

Next, to fix the notations, we identify the tangent spaces T Pi D rji with W l . It is well known that this 
identification yields normal coordinates centered at the points pi , namely 

x : B(p u 1) — > R" and y : B{p 2 , 1) — > E". 

We introduce now asymptotic cylindrical coordinates on the punctured ball B*(0, 1) = x (B(pi, 1) \ {pi}) 
setting t := loge — log |x| and 9 := x/\x\. In this way we have the diffeomorphism B*(0, 1) ~ (loge, +00) x 
S n_1 . Analogously, wc consider the diffeomorphism y (B(p2, 1) \ {pi}) = B*(0, 1) ~ (—00, — loge) xS n_1 , 
this time setting t := — loge + log \y\ and 9 := y/\y\. 

In order to define the differential structure of M £ , we excise a geodesic ball B(pi,e) from D m , obtaining 
an annular region A(pi,l,e) := B(pi,l) \B(pi,e), i = 1,2. The asymptotic cylindrical coordinates 
introduced above can be used to define a natural coordinate system on the neck region 

(t,9) : [%,l,e)U%,l, £ )]/ ► (loge, -loge) xS"-^:^, 



(> 



where ~ denotes the relation of equivalence which identifies the boundaries of B(p\,e) and B(p 2 ,e), 
namely 

gi~g 2 x/\x\(qi) =y/\y\(q 2 ) and \x\(qi) = e = \y\(q 2 ) . 

Clearly, in this coordinates, the two identified boundaries correspond now to the set {0} x S™ -1 . To 
complete the definition of the differential structure of M £ it is sufficient to consider the old coordinate 
charts on D rji \ B(pi, 1), i = 1,2. 

We are now ready to define on M £ the approximate solution metric g £ . First of all, we define g e to be 
equal to the gi on D ni \ B(pi, 1), i = 1, 2. To define g £ in the neck region, we start by observing that the 
choice of the normal coordinate system allows us to expand the two metric g\ and g 2 around pi and p 2 
respectively as 

.91 = [S a p + 0(\x\ 2 )] dx a ® dafi and g 2 = [6 af3 + 0(\y\ 2 )] dy a ® dy p . 

Recalling that the metrics g\ and 172 are locally conformally flat and using the (t, #)-coordinatcs introduced 
above, we can write 

gi = uf^(l + a)(dt 2 + d9 2 ) 1 with := e^e" 2 ^' 

92 = "2^(1 + c 2 )(dt 2 +d6 2 ), with u 2 {t) := e^e^* 
Now, we fix as background metric on M £ the following 

on D, H \ B( Pi , 1) 
on A(p lj l,e)l_lA(pa J l,e)]/~ 

where 

c := r/ci + (1 - r/) c 2 , 

with 77 a smooth and non decreasing cut off function such that 77 : (loge, — loge) — > [0, 1] and identically 
equal to 1 in (loge, —1] and in [1, — loge). Subsequentely, we consider another non increasing smooth 
function \ '■ (loge, — loge) — > [0, 1] which is identically equal to 1 in (loge, — loge — 1] and which satisfies 
lim t _j._i og£ x — 0. Using these cut-off functions, we can now define a new conformal factor 




1 onD m \B(pi,l) 
X(t) «i + x(~t) u 2 on A(p u 1, e) U A(p 2 , 1, e)]/ 



(4.1) u £ := 
Finally we define on M e the family of approximate solution metrics g 6l by setting 

4fc 

(4.2) g £ := uf^g . 

To conclude this section, we observe that with this definition we immediately have that for every m € N 
the approximate solution metrics converge to gi on the compact subset of D ni \ {pi} with respect to the 
C m -topology when the parameter e tends to 0, for i = 1,2. For these reasons, we expect that the size of 
the term Af(u £ ,g), which represents the fail of u £ from being an exact solution, will become smaller and 
smaller when e — > 0. Finally, we notice that adapting the proof of [4, Lemma 3.2] it is straightforward to 
show that for e sufficiently small g £ lies in ri,. 

5 Analysis of the linearized operator about the Delaunay-type 
metrics 

In this section we discuss some boundary value problems for the linearized operator introduced in (2.3) 
about Delaunay-type metrics. This local analysis, will find its application in Section 6. 
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We start by recovering from [17] the expression for the linearized operator about a Delaunay-type metric 

4fe/(n-2fe) w , 

9d,ti = v D ' v y g cy i. We set 

(5-1) h := vl„ - (^) 2 vl n and F := vf^/{H + «ff ) 

and we recall that that h(t) > for any t £ R (see [7] and [17]). With these definitions at hand, we can 
state the following 

Lemma 5.1 (Linearization about the Delaunay-type metrics). The linearized operator about the Dealunay- 
type solution Vn,r\ is given by 

(5.2) L {vD,n , 9cyi) W>] = - C n . k v D;V h~^~ { <9 t 2 + a v Ag - p n } [h^ w ] ; 

where Ag is the Laplace- Beltrami operator for standard round metric on gg™- 1 an d the coefficients a v and 
p v are given by 

( 5 ' 3 ) a n := k{n-l) + k[n-l) F ' 

(KA\ „ ._ ( n-2k \ 2 , n(nfc + n-2fc)(fc-l) p n 2 (fc 2 - 1) F 2 

w(2fcw-w + 2fc) „,;^k , n 2 fe(fc-l) „,«=fe 
4fe u £>,ij ^ H ^ 

and i/ie constant C n ^ is defined by C„ & := (tip) ^ 7 ' or rotational convenience we also define 

the conjugate linearized operator by 

(5.5) := c\ 2 + a,, A g - p n . 

Moreover, we have that there exists a positive constant c = c{n, k) > such that for every j > n + 1 and 
every admissible Delaunay parameter rj 

(5.6) a, ; Aj + p,, > c , 

where the positive real numbers \j, j £ N, denote the eigenvalues (counted with multiplicity) of Ag, i.e., 
-A e <pj = Xj <j>j. 

As a consequence of the last inequality, we will obtain the coercivity of the conjugate linearized operator 
Cr/ along the high frequencies (i.e., for j > n + 1). 



5.1 Jacobi fields 

Using the conformal cquivariance of the equation we introduce new families of solutions which are vari- 
ations of the standard Delaunay solution with neck-size parameter rj. The infinitesimal generators of 
these variations will provide us with natural elements sitting in the kernel of the linearized operator 
L(i>_d,t;, 9cyi) about UD,r/, namely the Jacobi fields. 

The first remark is that since the equation (3.2) is autonomous, then the solutions are translation invariant 
(with respect to the t variable). In particular, for r > 0, the functions Vjj „ T (t), defined by 

(5.7) VD,r,At) : = V D , v (t + l0g(T + l)) , 

are still solution to (3.2). To find other possible families of solutions it is convenient to use the conformal 
equivariance of equation (3.3). First notice that, writing t = — log |x| and 8 = x/\x\, with \ {0}, 

the cylindrical metric and the Euclidean one are related by g cy i = \x\^ 2 g^ on R™ \ {0}. As a consequence 
of (2.2) we get 

(5.8) M{v,g cy i) = \x\ n N{\x\~^v, m ~) . 
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Hence, if v(t, 8) solves Af(v, g C yi) = on R x §" 1 , then the function u(x), defined on R" \ {0} by 

(5.9) u(x) := \x\~ :L ^ h ~ v(— log \x\ , x/\x\) , 

is a solution to A/"(it, <7R") = on R" \ {0}. In particular, the Delaunay solutions VD. v (t) defined on the 
cylinder correspond to the radial solutions of the latter equation ud^(\x\) := | a; |-( n - 2fe )/ 2 ' c VD r) (— \ g | x |) 
with a pole in the origin. Since the equation satisfied by U£, :Tj is clearly translation invariant (due to the 
fact that the background metric is <?r>i), we have that, for b £ R™, the n-parameter family of functions 
U£>,n,b( x )i defined by 

(5.10) UD, v ,b( x ) := ud, v (\x - b\) = \x-b\ ^~vd, v (- log \x - b\) , 

still satisfies N(ud^^ ,gw) = 0. These functions present a singularity at b £ R™ and they are radial 
with respect to b £ R™. These new solutions UD,rj,b{x) defined on R n \ {b} correspond via (5.9) to the 
solutions v Djhb (t, 6) of the equation (3.1) defined on R x §' l_1 \ {(- log |6| , 6/|&|)} by 



(5.11) 



VD, n ,b(t,0) 



\6 - be^—^r-VD^t - log \6 - be 1 



The last family of solutions comes in the following way. First observe that the function t i— > VD,r/{t) '■= 
VD,r](—t) is still a solution to (3.1) on R x § n_1 . This corresponds to the fact that on R\ {0} the Kelvin 
transform of md.,, namely the function 

n-2k ., |OU , n — 2k 



(5.12) 



UD,r,(\x\) 



V D ,r,{- log \x/\x\* 



satisfies Af (md,^ <?r™) = 0. Now we translate UD tV by a vector a £ R", obtaining an n-parameter family 
of functions UD,rj,a(x) := U£,_ v {\x — a\) = \x — a|~(™~ 2fc )/ 2fc V£). r) (log \x — a|) and finally we take the Kelvin 
transforms of the Uo,ri,a s obtaining, for a £ R™ the new family of solutions on R" \ {0} 



(5.13) 



UD,n,a(x) 



X 



a\x\ 



v D 



2 log \x\ + log la: 



a\x\ 



These solutions are no longer radial and present a singularity at the origin. For a £ 
on R x S" -1 \ {(log |o| , a/\a\ )} to the solutions 

- tI — 



they correspond 



(5.14) 



ae 



\~~ ss ~v D , fl (t + log |6» - ae-*|) 



In the remaining part of this section we will use all these families of solutions to define some special 
elements in the kernel of the linearized operator around the Delaunay solutions ■ First of all, we 
recall that if A i— > i>.d, j; ,a is a variation of vd i?) such that for every admissible value of the parameter A 

AT (vD, v ,\,9cyi) = and v D ,r,,o(t) = v D , v (t) , 
then it is straightforward to see that 







d_ 



A/" (UD,,j,A j 9cyl) = M V D^ , g C yl) 



X=0 



d_ 



VD, ; 



A=0 



where TL(vd.ti ;9cyl) represents the linearized operator around the Delaunay solution vd.t}- The functions 
^a|> = o v D,r],\ ar e the so called Jacobi fields and they clearly belong to the kernel of h(i)D iTI , g cy i)- Applying 
this reasoning to the family of solutions a H> WD )TJ + a and r i-> Wd )T7)T , it is natural to define the quantities 



(5.15) *»•-(*) 



d_ 

da 



v D,ii+a{t) 



and 



n 



(0 



d_ 



VD,r,,r(t) = *D,jj(*) 



In analogy with that, we use the other two families b i— > vu,ri,b and a n- VD,n,a to define, for j = 1, 
the Jacobi fields 

<9^ 



(5.16) 
(5.17) 



VD, n ,b(t,6) 



b=0 



~ n-2k 
. 2k 



VD,<n(t) + i'D, n {t) 



M0) , 



9aJ 



VD,r),a(t,6) 



\ n-2k 
V 2k 



VD,ri(t) ~ V D , v (t) 



a=0 



where the 0j's are the n eigenfunction of the Laplacian on S n with eigenvalue n— 1, namely — A 
(n- 1)0 3 , for j = l,...,n. 
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5.2 A linear problem on the cylinder E x § n 1 

In this subsection wc want to study the problem 

L(v D!V ,g cy i)[w] = f inKxS"- 1 . 

Following [16] we observe that the natural functional setting for this problem is given by weighted Holder 
or Sobolev spaces. Both choices are essentially equivalent. However, in our argument we will use Holder 
spaces. For a fixed weight parameter 6 £ K and m € N we define the space 

Cf{D v ) := {ueC m (D v ) : \\u\\c r < +^} , 

where the weighted norm is defined by 

IMIcyCAO : = SU P E£i(«*h*)~ f |V J 'u|(t,0). 

KxS"" 1 

We point out that | • | and V are respectively the norm and the Levi-Civita connection of the cylindrical 
metric g cy i- In the same way we define for 5 € K, m G N and j3 € (0, 1) the weighted Holder seminorm by 

(5.18) [u] c m.p (D ) := sup (coshi)"" 5 [u] C m,3 ((t _ 1 t+1)xS „-i ) . 
The weighted Holder spaces are then given by 

(5.19) CT'^D,) := {ueC m ^(D v ) : Hl^,, := \\u\\c r + [n} c ^ < +oo } . 
Following the analysis in [18], one immediately find that 

Hv D , v ,gcyi) ■■ c 2 /(d v ) — ► Cj l/S (D,) 

is Fredholm, provided 5 ^ I v , where := {±Sj,ri ■ j G N} is the set of the indicial roots of the operator 
L(vd,t;, 9cyi) at both +oo and — oo. In general the indicial roots (for a precise definition see [18]) do 
depend on the neck-size parameter r], but here it follows from the explicit knowledge of the Jacobi fields 
that Jo, 77 and Si^ n are independent of 77. In particular the indicial root 5o,n = is related to the Jacobi 
fields and ^ / ° ,+ , which are respectively linearly growing and bounded in t, whereas the indicial 

root 5i^ n = 1 has multiplicity n and is related to the Jacobi fields ^F^' - and j = 1, . . . ,n, which 

are respectively exponentially growing with rate e* and exponentially decreasing with rate e . We also 
point out that as a consequence of the inequality (5.6) it can be deduced that for every admissible value 
of the Dclaunay parameter 77, the indicial root 82, rj verifies the inequality 

(5-20) S(n,k) := + (^f < S 2 ,, . 

We are now in the position to prove the following 
Lemma 5.2. Let 1 < 5, then the operator 

HvD, v ,g cv i) ■■ C 2 J(D V ) — ► C°J(D V ) 

is injective. 

Proof. Since the functions which are involved in the conjugation (5.5) are bounded and positive, it is not 
restrictive to prove the result for the conjugate operator 

C v : C 2 J(D V ) — ► C°J(D V ) . 

Performing a standard separation of variables and projecting the equation along the eigenfunctions of Ag, 
we note that for the low frequencies j = 0, . . . , n the space of the general solutions to the homogeneous 
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equation is spanned by the (conjugate) Jacobi fields $^ := hS k 1 '/ 2 5 , ^ ,± . On the other hand, it is easy 
to check that, for 5 > 1, no one of these functions belongs to the weighted space C^(D V ). Thus, it 
is sufficient to test the injectivity only for the high frequencies, j > n + 1. Hence, suppose to have a 
function <I> such that 

£„$ = and $M) = Ej> n+ i * J '(*)^W • 
Since $ G C^f (A,), wc have that 

1*1 (M) < C- (cosht)" 15 

for some fixed C > 0. On the other hand, the maximum principle (which holds when C v acts on the high 
frequencies, see [17]) gives 

|$|(T,0) < C- (cosh TV 5 
for every T £ R. Letting T — > +oo we deduce that $ = and the proof is complete. □ 

Using the fact that C v is formally selfadjoint, it is standard to deduce (see [16]) that 

h(v D , v ,g cyl ) : C 2 /(D V ) — > C°/(D n ) 

is surjective for 5 > I, 6 I n . Following [16] we are going to improve these first issue by showing that 
the surjectivity can be obtained on a smaller space. To do that it is convenient to set 

W(D V . R ) := span{xR^' ± : i = 0,...,n} , 

where \R is a non decreasing smooth cut-off function which is identically equal to 1 for t > R and which 
vanish for t < R — 1. A simple adaptation of the ODE argument used in [16, Proposition 2.7], gives us 
the following 

Lemma 5.3. Let 1 < 5 < 6(n,k), then the operator 

Hv D , v ,g cyl ) : C 2 J(D V ) © W(D, hR ) — ► C°J(D V ) 

is surjective. 

5.3 A linear Dirichlet problem on the half cylinder R + x § n_1 
In this subsection we study the Dirichlet problem 

jh{v D:V ,g cy i)[w} = f in (R , +oo) x S"" 1 , 
1 w - on {R} x S"" 1 , 



(5.21) 



for which will prove a well posedness result in the next Proposition 5.6. As it will be apparent from 
the proof, this result heavily relies on a proper choice of the value of R. Loosely speaking, the correct 
choice of R has to compensate the lack of maximum principle for the linear operator L(wo. ?? , g cy i). The 
same kind of problem will show up also in the next subsection 5.4. For future convenience we set 
Dr/,R '■= (R • +oo) x §™ . Again, since h and vd,t/ arc bounded and periodic we study without loss of 
generality the conjugate problem 

(5.22) (£r,z = V in A,,*, 

z = on dDr, tR , 

where z = h t ^ k ^ 1 ^ 2 w and y = ~Cn\ v ~DrJ ri ~ i ' k ~ 1 ^l' 1 f ■ Wc consider now the usual eigenfunction decompo- 
sition 

oo oo 

y(t,B) = Xy'(t)&(0) and z(t,6) = £V(t) , 

3=0 3=1 
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where the (fij's indicate the eigcnfunctions of the Laplacc-Bcltrami operator on (S n ~ ,gs"—i ) which 
satisfy the identities — Ae <\>j = Xj 4>j , with j £ N. We also recall that the spectrum of Ag is given by 
{ m (n — 2 + m) : m £ N } and that in particular the first nonzero eigenvalue is n — 1, with multiplicity 
n. 

In the spirit of [15], it is convenient to treat separately the high frequencies, i.e., j > n + 1, and the low 
frequencies, namely j = 0,...,n. Basically, this distinction is motivated by the fact that, depending on 
the size of the Xj , the quantity a v Xj + p v presents a change of sign and this has a clear influence on the 
analytical properties of our operators. 

High frequencies: j > n + 1. We consider the projection of z and y along the high frequencies 

00 00 
y(t,6) := J2 V* and Ht,6) ■= E ^(*)&(*). 

j— n+1 j—n+1 

and for T > R we consider the projected and truncated linear Dirichlet problem 



(5.23) 



C v z = y in D^ R , 
z = on dD^ R , 



where := (R,T) x 8" 1 . In the high frequencies regime the linear problem (5.23) has a clear 

variational structure. Indeed it is easy to see that critical points of the Eulcr-Lagrange functional 

(5.24) E T {z) := f [ (\d t z\ 2 + a v \V g z\ 2 e + p v z 2 + 2yz )dtd6 

Jr Js"- 1 

are weak solutions of (5.23) (here d9 represent the volume element of the round metric gs™- 1 on the 
(n — l)-dimensional sphere). On the other hand, since j > n + 1, we have by [17, Lemma 5.3] that 
a v Xj + > in D n R. This implies that the functional Et is coercive on 

[ Hi (Dl R ) f.Lflj (Dl R ) J u(-,e)^j(6) d6 - , j = 0, . . . , n 

hence it is bounded from below. Furthermore it is easy to check that the functional Et is weakly lower- 
semicontinuous on [Hq [LiE R )\ . Thus, using the direct method of calculus of variations, we infer the 
existence of a minimizer zt of Et, which provides a (weak) solution of (5.23). The standard elliptic 
theory yields the expected regularity issues for zt in terms of the regularity of y. 

Moreover, as a particular case of [17, Proposition 6.4], in the high frequencies regime, there holds the 
following 

Lemma 5.4. Let \ S\ < 6(n,k), then there exists a positive constant C = C{n, k,S) > such that if 
zt G C^g (D^ R ) and y € Cjg (D'E R ) verify (5.23), then we have 

(5.25) || z T \\c^ (Dl R ) ^ C WyWc^(Dl R ) . 
for every T > R. 

Using the fact that the estimate is independent of T, it is easy to obtain a solution z to (5.22) by letting 
T — > +00. Moreover it is clear that z verifies the estimate 

(5-26) \\ 2 \\c^(D n , R ) < C WvWdhUD*,*) ' 

with the same constant C as in Proposition 5.4. 
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Low frequencies: j = 0, . . . , n. Here we start by considering the projection of our original problem 
(5.22) along the eigcnfunction cj>o, obtaining 



(5.27) 



Z v ° z° = y° in ( R , +oo ) , 
z°(R) = . 



where £® := &f — p n . As it is evident, in this case the potential has a wrong sign, thus we are forced 
to use a different approach in order to provide existence. We suppose that the right hand side is at least 
continuous and we extend it to the whole R (with a small abuse of notations, we still denote this extension 
by y°). Next, following [15], we consider, for any T > R, the auxiliary backward Cauchy problem 

(5.28) 



z 


= y° 


in ( — oo , T 


< z(T) 


= o , 




. HT) 


= . 





Using the Cauchy-Lipschitz Theorem, we infer the existence of a unique solution Zj. to (5.28). As we are 
going to show, the weighted norms of these solutions admit a bound which is uniform in T . This will 
allow us to produce a solution to the problem (5.27) with the wrong boundary data, just by taking the 
limit of the z^s for T —> +oo. As a final step we will correct these boundary data by adding a suitable 
multiple of the (conjugated) Jacobi field $° ,+ := ft,( fe_1 )/ 2 = frC^ 1 )/ 2 vd, v , which lies by definition 
in the kernel of L ° . 

Lemma 5.5. Let < <5 and R = rhT ri +r with rh £ N and f e 1 sufficiently small. Then, there exists 
a positive constant C = C(n,k,5) > such that if z^ £ C^g (R,T) and y° € C^g (R,T) verify (5.28), 
then we have 

(5.29) H4ll c ^( fi) T) < C \\y°\\c°j§«R,T) . 

for every T > R. 

Proof. We only prove the weighted C°-estimate, since the weighted C 2, ^-estimate will follow by standard 
scaling arguments. We want to establish the T-uniform bound 

II z t\\c° s (r,t) < C II V \\c°_ s ((r,t) ■ 

We argue by contradiction. If the statement does not hold, then it is possible to find a sequence of triples 
(T^ZT^Vi) such that 

• £°4. = Vi in (R, T i) and ^.(Ti) = = z\{Ti) for every i G N , 

• II z Ti \\c° 6 (R,Ti) = 1 for ever y « e N , 

• II Vi \\c° s {R,Ti) — >0 as i-» +oo . 

From the second point we infer the existence of a point U £ (i?,T, ) such that 

sup e st \z^\(t) = e su \z^\(U) = 1 . 
te(R,T t ) 

In fact on the half cylinder the weighting function cosht can be replaced by e* in the definition of the 
weighted norms. This yields an equivalent norm and simplify the computations of this subsection. It is 
now convenient to set, for every i £ N, 

Zi{t) := e su z^(t + t t ) and Vi (t) := e St ' y? (t + U) , 

where the point t varies now in ( — logi? — ti , T,; — ti), for every i £ N. From these definitions it follows 
that 
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£° gi = Vi in (R-ti,Ti- U) and z,(Ti - U) = = i" (T< - t<) for every i£N, 
sw Pte(R-n .Ti-u) e5t \ z i\ (t) = N(0) = 1 for every ieN, 



• sup te( _ fi _ t . iT ._ ti) e 5 * 1 2/i| (t) — ► asi^+oo. 

We are now ready to let i — > +oo and study the different limit situations in order to get a contradiction. 
As a first step we remark that (up to a subsequence) the intervals (R — ti, Ti — 1{) converge to an interval 
(/3 _ , (3 + ) which is nonempty. In fact, since R — ti < and Ti — ti > 0, we have immediately that 
f3~ 6 R~ U {— oo} and f3 + £ R + U {+oo}. Moreover, we claim that (3 + is strictly positive. In fact, if it 
would not be the case, then we would have that up to a subsequence Ti — U — > 0. Since | Zi\ (0) = 1 
and Zi {Ti — ti) = = z% (Ti — ti) for every i g N, the quantities | dtZi |'s must explode in the intervals 
( Ti — ti — 1 , Ti — U ), as i — > +oo. On the other hand, from the hypothesis on z% and j/j it follows easily 
that 

\d? Zi \(t) < Ce- 5 ^-^, 

on the interval (Ti — ti — 1 , Tj — ti). Since we are supposing that Ti — ti — > 0, this inequality tells us 
that the second derivatives of z% are uniformly bounded as i — > +oo. The fact that dtZi (Ti — ti) = for 
every i £ N implies that the first derivatives of Zi also admit a uniform bound on ( Ti — tj — 1 , Ti — ti ) 
as i — > +oo, which is a contradiction. Hence we have that (/3~, is always nonempty. 

The equation satisfied by the z^s implies that there exists a function uJqo such that Wi — > -u^ in 
C^ oc (f3~ , /3 + ). In particular, the function z^ verifies the homogeneous equation 

(5.30) £°Zoo = in (/T,/3+), 

in the sense of distributions. As a consequence z^ can be written as a linear combination of the Jacobi 
fields and $°' + , namely, there exists A, B e R such that 

Moreover, the hypothesis on | z.i \ (0) implies at once that | z x \ (0) = 1. Thus z^ is non trivial. When 
(3 + < +oo, then the Cauchy data for the limit problem are given by z^ (/?+) = = ioo ((3 + ), thus 
Zoo = and we have a contradiction. If /3 + = +oo, then the decay prescription | Zoo\(t) < e~ St with 
<5 > implies that both the constants A and B must be zero, contradicting the non triviality of z^. □ 

Since the estimate (5.29) is independent of the parameter T > R, we let T — > +oo and we obtain a 
function z° which verifies the identity 



C° z° = y° in (R,+oo) 



together with the T-uniform estimate 



(5-31) P U lkf(fl,+oc) < C lly U llc»f(«, + oo) . 

where 5 > and C is the same constant as in Lemma 5.5. 

The next step amounts to correct the function z° to a solution of the problem (5.27). This will be done 
by adding an element in the kernel of L ® to the function w in order to fulfill the homogeneous boundary 
condition at t = R. Here we decide to choose the (conjugated) Jacobi field 3>°' + := h^^'^VD,^- We 
notice that this correction is no longer a function in C^'g(R , +oo), with 5 > 0, since it is just bounded 
at +oo. With these considerations, we are now ready to set 

(5-32) At) := z\t) - J££L . 

<P V (R) 
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It is now immediate to check that this yields a solution to (5.27). We point out that the definition of w° 
makes sense since R = mT n + r and thus 7^ 0. From the definition of z° it follows at once 

that its component along the Jacobi field $°' + is bounded by (C/$°' + )(i?) || y°\\co,p^ R +oa < ) with C as in 

Proposition 5.5, hence for S > the solution z° to (5.27) is unique in the space 

Clf (R,+oo) ® span {$"<+}. 
Now we are ready to treat the projection of (5.22) along the eigenfunction cf>j, with j = 1, . . . ,n 



(5.33) 



z j = y j in ( R , +00 ) , 
w j (R) = , 



where £i := <9 f 2 — Xja v — p v . Proceeding in the same manner as in the case J = we deduce that for 
6 > 1 there exists a unique solution z J to this problem in the space 

C*f{R, +00) © span{^+} 

which can be written as 

* J = *> + TT=5E 7±^~. 77^T^ + , J = l,...n . 



2k 



: -VD.r,{R) - VD,n(R) 



Note that this definition makes sense since, thanks to h(t) > for any t GM. (see (5.1)), there holds that 
[ n ?k k vd,t)(R) — vd.ti(R) ] 7^ 0. Moreover we have as in the previous case that z? verifies the estimate 

( 5 - 34 ) M j Hcif(H 1+ oo) ^ C Wy 3 Wc°:l{R,+oo) > 

where now 5 > 1 and C is a positive constant only depending on n, k and S and the component 
of z j along the (conjugate) Jacobi fields := /iO" 1 )/ 2 ^ is bounded by ( C / [ ^^Vd^R) - 

*D,u(-R)]) l|y J 'llc°f(R,+oo)- 

To summarize all the result of this subsection, we define the finite dimensional function space 

W + (D VtR ) := span{^+ : j = 0,...,n} 
and for a function u = Y^j=o a t + e W + (D, h fi) we simply set 

(5-35) \\ u \\w+(D v , R ) ■■= E"=ol a |l- 

We thus have proved the following 

Proposition 5.6. Let 1 < S < S(n, k) and R = rhT v + f as above, then for every f 6 C^g (D^^r) there 
exists a unique solution w G C 2 .'^ (D^^r) (BW + (D r)t fi) to the problem 

j h(v Dir ,,g C yi) [w] = f in (R , +00) x §" _1 , 
\ w = on {R} x S"- 1 . 

Moreover we have that there exists a positive constant C = C(n,k,S,T)) > such that 

IHI C =f(A,, H )ffiW+(A,, R ) := Mc'jg + IMIw+CA,.*) < C \\f\\c^(D 71 , n ) ■ 

The same analysis can be reproduced on a domain of the form D, h -R := (—00, —R) x §™ _1 , in order to 
solve the problem 



(5.36) 



M v D, n ,g cy i)[w} = f in (-00, -R) X S n 1 
10 = on {-#} x S"- 1 , 
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where R is a real number of the form R = rhT v + f , as in Proposition 5.6. In this situation we use the 
finite dimensional function space 

W~(D,,_ji) := span {^/- : j = 0,...,n} 

for the corrections along the low frequencies. Obviously, for a function u = Y^j=o a J ^r}~ ^ W + (D, h _n) 
the norm can be defined by 

II "II w- (£>„,_*) := Ej=ol a 7l ■ 
In analogy with the previous proposition, it is straightforward to obtain the following 

Proposition 5.7. Let 1 < 5 < 5{n,k) and R = mT„ + f as above, then for every f £ C^'f (D„ _r) 

f/iere exists a unique solution w £E (D^^r) W _ (Z)^ i _fl) £o i/ie problem (5.36). Moreover we have 
that there exists a positive constant C = C(n,k,5,rf) > such that 

IMIcif(A,,_ H )ew-(^,_ K ) := IMIcif (d„,_ r ) + IMIw-(A,,_ fl ) < C H/ll c °f(D,,_ R ) ■ 

To conclude we observe that as a consequence of the analysis of this subsection we can also solve the 
same problems with nonzero boundary condition, namely for every 1 < 5 < 5(n, k), every / e C°^g(D v ) 

and every boundary datum <E C 2,!3 (dD v .±u), there exists a unique function <E C^g (D Vt ±u) 
W ± (D v ,±r) verifying 

|LK„flc ! ;i)[w ± ] = / in A,,±i? , 
ui = v on dD rh ±R , 

together with the estimate 

H w± llcif(^, ±R )ffiW-±(B,, ±B ) ^ C [ H/ll C°f(A,, ±B ) + \\ v± \\c^{dD v , ±R )] , 
for some positive constant C = C(n, k, 8, rj, R) > 0. We thus have proved 
Proposition 5.8. Let 1 < 6 < 6 and R = mT n + f as above, then the operator 

Uv Dtr „9cyi) ® d ± : C«if (A,,±fl) © W ± (^„,± fi ) — > C°f (A,,±ii) x C 2 ^ {dD 1h±R ) 



is an isomorphism, where : w i— > w\dDr, ± R - 



5.4 A linear Dirichlet problem on a finite cylinder (— i2, i?) x § n 1 

In this subsection we study the Dirichlet problem 

[L(t; A „yH = / in (-i?,i?)xS"- 1 , 
(5.37) < 

[ w = on x S"- 1 , 

for which we are going to prove the following 



Proposition 5.9. Let R = mT, + r with rh £ N sufficiently large and f € R sufficiently small. Then, 
for any f e C°'P((-R, R) x S n_1 ) iftere exists a unique solution w G C 2 'P((-R,R) x S™" 1 ) £o (5.37). 
Moreover, there exists a positive constant C = C(n, k) such that 

(5-38) IMIc^u-fl^xS™- 1 ) < C ll/llc°.' 3 ((-i?.,i?,)xS"" 1 )- 
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Proof. Thanks to the compactness of the domain, we use the Fredholm alternative to prove the well 
posedness of (5.37). Thus, the existence and uniqueness of solutions to (5.37) follows from the fact that 
the homogeneous problem 

(Uv Difl ,g cyl )[w} = in (-R,R)x§ n - 1 , 
\ w = on {±R} x S™- 1 , 

only admits the trivial solution w = 0. This is equivalent to say that z = is the unique solution of the 
conjugate problem 



(5.39) 



£ v z = 
z = 



in (-R,R) x S"- 1 
on {±R} x S"- 1 , 



where z := h^- k ~ 1 ^ 2 w. To show that z = uniquely solves (5.39) we project the equation along the 
cigcnfunctions 4>j of the Laplace Beltrami operator on the sphere (§™ _1 , gg™- 1 )- Then, as we already 
did in the proof of Proposition 5.6, we treat separately the high frequencies (J > n + 1) and the low 
frequencies (j = 0, . . . , n). Thus, we decompose z as 



(5.40) 



3=0 j=n+l 



Now, since in the high frequencies regime the (weak) solutions to (5.39) can be obtained as critical points 
of the coercive and weakly lower semicontinuous energy defined in (5.24), it turns out that the second 
sum in (5.40) is identically equal to zero. 



n and we note that 
•j(0), where the 



To obtain the same result for the low frequencies, we start with the case j = 1, . . 
the Fourier coefficients z J can be written as a linear combination of <£>^ + 4>j(ff) and 
conjugated Jacobi fields are defined as := h^ -1 '' 2 ^^. Hence, there exist real numbers Aj 

and Bj such that 

**(t) = h^-^it) {A 3 [^v D . v {t) + « A ,(t)]e* + B 3 -[2=^t; Di)J (i) - i D ,,(t)] e" 4 } . 

We are going to show that the homogeneous Dirichlct boundary conditions together with our choice of 
R imply that all the z' ] must vanish for every j = 1, . . . ,n. From the null boundary condition we deduce 
that 




VD,t)(R) 



tanh (R) 



n-2k 

2k " v D , n (R) 



2k 



tanh (R) 



VD,T){R) 



= 0, 
= 0, 



where we have used the fact that «£>,«(•) is an even function. Now, using the fact that R is of the form 
R = mTr) +f, it is sufficient to choose fa large enough and f small enough in order to insure that 

j_2fe , VD,v( R ) 



2k VD, V (R) 

(R) 



2k 



tanh (R) 



tanh(i?) ^ 



0. 
0. 



Hence, A,- 



VD,r,(R) 

= Bj, for every j = 1, . . . , n. 
It remains to prove that also z° = 0. Now, z° has this form 

z °(t) =ft(*- 1 )/ 3 (t) {A)*°'-(t)+Bo*S' + W} 



17 



for suitable Aq and Bq in R. From the homogeneous boundary conditions combined with the fact that 
is odd and is even, we deduce at once that if R is chosen as in the statement both A and B 
must vanish. 

Hence, the problem (5.37) is well posed. On the other hand the a priori estimate (5.38) is a direct 
consequence of the standard elliptic regularity theory. □ 



6 Global linear analysis 

The aim of this section is to provide existence, uniqueness and a priori estimates for solutions to the 
linear problem 

(6.1) h(u B ,g)[w] = f in M £ . 

For sake of simplicity we just consider the case where M e is the connected sum of two Delaunay 
type solution M E = D m %D m . All the arguments can be trivially adapted to the general case M e = 

Ajitte ■ ■ 4eD VN - 

To introduce the correct functional framework for the global linear analysis on M £ we define 

(6.2) ll"llc^"(M«) := IMIc«-''(A, 1 \B(p 1 ,l)) + ll U !lc7^(A J2 \B(P2,l)) + NIc^JV.) . 

where the first two norms are defined as in Subsection 5.2 and the norm in the neck region N E is defined 
by 

(6.3) \\u\\ c ^ (N) := sup £7 =1 ( £ coshi) 7 |V^|(^) 



+ sup (ecoshi) 7 [w] c . m ,^ ((t _ lit+1)xS „-i ) . 

t e (logs,- loge) 

Note that, again, | • | and V are respectively the norm and the Levi-Civita connection of the cylindrical 
metric g cy i whereas [ • ] stands for the usual Holder seminorm. As a consequence, we introduce the 
following weighted Holder space 

C?,f(M e ) := {ueCr o f(M e ) : \\u\\ c ^ {Me) < +^}. 

6.1 Global and uniform a priori estimates on M £ 

We recover from [4] the following removable singularities result 

Lemma 6.1 (Removable singularities). Let g = (1 + b) ™- 2fc <7r™ be a conformally flat metric defined on 
a geodesic ball B(p, 1) verifying the equation 

a k (g^A g ) = 2- k (l) . 

Suppose w satisfies in the sense of distributions 

L(l + b >m n)[w] = on B*(p,l) 

with \w(q)\ < C\distg{q,p)\~ tl for any q £ B*(p, 1) for some positive constant C > and for some weight 
parameter < [i < n — 2. Then w is a bounded smooth function on B(p, 1) and satisfies the equation 
above on the entire ball. 

We are now in the position to prove the following e- uniform a priori estimate for solutions to (6.1). 
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Proposition 6.2. Suppose that 1 < 5 < 8 (n, k), < 7 < (n — 2k) /k and let w G C 2 '^ „_ 2 k (-^e) 

— °>7 5TT" 

/ G C°'f 7 _( n _2fc)(-^e) ^ e ^ zt ' functions satisfying 

h(u e ,g) [to] = / in M e . 

Then there exist C = C(n, k, /3, 7, <5) > and £0 = £o(n, fej 7j 5) sucft that, for every e G (0, £o], we have 
(6.4) \\w\y,, < C ||/|| c o,. 

Proof. Before starting the proof, we advise the reader that we will prove (6.4) using a different, albeit 
equivalent, norm. With a little abuse of notation, we introduce the norm 

\\ U \\c™£(M e ) ■■= maX {\\ U Wc™f(D vl \B( Pu l)y W U \\c^(D V2 \B( P2 S)V IHc^ (N e )} > 

which is clearly equivalent to (6.2). We will just provide the uniform weighted C°— bound 

Hie- 2k < c ||/|| c o , 

n Ik 

since the uniform weighted C 2 '^-bound will follows easily from standard scaling argument, sec [18]. To 
prove the above estimate we argue by contradiction. Suppose that there exists a sequence (£j, Wi, fi), 
i G N, such that 

• £j — > 0, as i ^ +00, 

• IHIc° si, i G N, 

• ll/i||c° , — asi^+oo 
and 

L(u £i , g) [wi] = ft in M e . 

Now, up to a not relabelled subsequence of i, one has to face with the following two cases: 

!• \\ w t\\c^ 5 (D^\B( Pj ,i)) = 1, for any i G N and for j = 1 or 2. 

2. ||wj|| c o(jv e ) = 1 for any i G N. 

The second case has been fully analysed in [4, Proposition 4.4, case 2. and case 3.] to which we refer for 
the details. Concerning the first one, we note that there is no loss of generality in restricting the analysis 
only to D Vl (recall that on D m we use coordinates r% and 6, according to Section 4). Secondly, it is 
natural to split the case 1 into two subcases. The first subcase appears when there exists M > and a 
subsequence of points qi = (rj,6>j)'s such that r t G [—A/, A I] x § n_1 and (cosh r i) 5 \wi\(ri, 9 i) > 1/2. The 
second subcase is when the points qi's leave every compact set of D Vl . However, it follows from (5.37) 
that this second subcase can always be reduced to the first one. The assumption qi G [—M,M] x S™ _1 
implies that, up to a not relabelled subsequence, there holds that qi — > q^, iUj — > lUoo in Cf oc (D Vl \ {pi}) 
and fi — > in Cf oc (D m \ {pi})- Thus, the function Woo is nontrivial, since |iUoo|(ooo) > (coshM)~ 5 /2, 
and solves in the sense of distributions the limit problem 

L(l,.gi) [u^ 1 Woo ] = in D m \{ Pl } with || u> o||c^(A, 1 \B(pi,i)) = 1> 

where u\ is the function used in the construction of the approximate solutions (see Section 4) , which we 
have set to be equal to 1 in D m \ B(pi, 1). If we show that the limit problem is solved by Woo on 
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the whole D Vl , then, using Lemma 5.2 (which is in force thanks to the fact that 1 < 6), we will reach the 
desired contradiction. To remove the singularity at p\, we observe that on B(pi, 1) we can always write 

,9i = (1 + 6i)"- 2fc 5R«, 

with bi(q) — 0(\dist gi (q,pi)\ 2 ). Hence, thanks to the conformal equivariance property (2.4), we have 
that the limit equation can be rewritten as 

L ((1 + h), 5R .) [(1 + 6iKV»] = on fl(pi,l)\{pi} . 
Recalling that g\ solves the cr^-Yamabe equation and that 

|(1 + b 1 )ui 1 w ao \(q) < C \dist gi (q, Pl )\-^ , 
we can apply Lemma 6.1 to obtain that u^ 1 w oa extends through pi to a nontrivial smooth solution of 

L(l,3i) [u^ 1 w 00 ] = in D m . 
This completes the proof. □ 
As a consequence of the a priori estimates, we obtain the following 

Corollary 6.3. Suppose that 1 < 8 < 5 (n, k) and < 7 < (n — 2k) /k, then there exists a positive real 
number £o(n, k, 7, 6) > such that, for every e € (0, eq] the operator 



L(«e,5) : C 2 ^_^(M e ) C°!^_ {n _ 2k) (M e ) 



is infective. 



The next task is to provide surjectivity for L(u £ , <?). Unfortunately, with this choice of the weight param- 
eter, which will turn out to be suitable for the nonliner analysis, the surjectivity cannot be recovered. 
The duality theory would suggest the use of the spaces with weight parameter S insted of —S, but as it 
is remarked in [16], these spaces are definitely too large. In particular they contain functions which may 
grow too fast at ±00 and even worst which are not integrable, in the sense which is made precise below. 
To overcome this difficulty, one can take advantage of the Fredholm character of the operator (which 
is actually the case, when 5 is not an indicial root) by considering a finite dimensional extension of the 
domain, the so called deficiency space. Of course, there are several different options for the choice of such 
a space (for example a different approach is contained in [16]). We start by introducing the following 
spaces 

W(D m >Rl ) := span { X R[ %f ■ j = 0, . . . , n } and W{D Vu - Ry ) := span { x~r[ %f ■ j = 0, . . . , 
W(D V2iR2 ):= span { XR ' 2 %t : j = 0,...,n} and W(D V2 _h 2 ) := span { X -R' 2 ^ : j = 0, . . . , 

W+ (D m :+Rl ) := span { X +i?; : j = 0, . . . , n } and W+ (D V2 :+R2 ) := span { X +r> 2 ■ j = 0, . . . 

W-(D Vu _ Rl ) := span{x- fli : j = 0,...,n} and W{D V2> _ R2 ) := span{x-^ : j = 0, . . . 

where Xr[ is a non decreasing smooth cut-off function defined on D m which is identically equal to 1 for 
t > R[ and which vanish for t < R[ — 1 (with — 1 > R\. The cut-off function X-R[ i s defined by 
the relationship X-R' 1 { r ) = Xi?i( — r )i an( i X-RJ, an( i X-i?a are defined on D m in an analogous fashion. 
Moreover, the parameters R[ and R' 2 are chosen in such a way that R[ — 1 > i?i and i?2 — 1 > R2, in 
order to apply the analysis of the previous section. 

We observe that all the functions in these spaces are integrable at ±00 in the sense that they are linear 
combinations of infinitesimal generators of families of conformal variations of the original Delaunay type 
solutions a nd vd,t/ 2 - This fact has an important geometrical meaning which will be made clear later 
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and which will be used in the nonlinear framework to justify the choice of a correction w with components 
lying in these spaces. To continue, we define the full deficiency space W(M £ ) by 

W(M e ) := W(D VltRl ) ® W(D m _ Rl ) ® W(D V2tR2 ) ® W(D m - R2 ) . 

We notice that since all of these spaces are finite dimensional, we can choose any norm on them. To be 
definite we always consider the norm given by the sum of the absolute value of the components. 

The importance of these spaces for the linear analysis is clarified by the following Proposition, which can 
be deduced combining Corollary 6.3 with [18, Therocm 12.4.2], 

Proposition 6.4. Suppose that 1 < d < S (n, k) and < 7 < (n — 2k) /k, then there exists a positive real 
number £q(u, k, 7, 5) > such that, for every e £ (0, £q] the operator 

h(u £ ,g) : C 2 J^_^(M e )®W(M e ) — ► C^ 7 _ (n _ 2fc) (M e ) 

is surjective and 

dim Ker L (u e , g) = \ dimW(M £ ) = 4(n+l). 
Our deficiency space is defined by 

(6.5) W(M e ) := W + (D VuRl )®W-(D m> . Rl )®W + (D V2>R2 )®W-(D ri2> - R2 ) . 

Incidentally we note that in [16], due to the use of a different functional framework, the deficiency space 
is chosen to be 

W(M £ ) := W(D VuRl )(BW(D V2 . R . 2 ) . 
As expected, dim W(M e ) = 4(n+l) = dim W(Af e ). 

The remaining part of this section will be devoted to prove the core of our linear analysis, namely the 
following 

Proposition 6.5. Suppose that 1 < 6 < S (n, k) and < 7 < (n — 2k) /k, then there exists a positive real 
number eo(n, k, 7, S) > such that, for every e £ (0, £0] the operator 

Uu s ,g) : C 2 ^_ 2 _ 2k (M s )®W(M s ) — ► C°J^_ (n _ m (M s ) 
is an isomorphism. Moreover the following e-uniform a priori estimate is satisfied 

(6.6) HHIc 2 '* 3 (M c ) + IMIw(M £ ) < C WL(u s ,g)[w]\\ c o,e (M) , 

where the positive constant C only depends on n, k, /3, 7 and 6. 

Proof. First of all, Proposition 6.4 furnishes for any / £ C°'f 7 _(„_ 2fe )(-M e ) the existence of a solution u 
to (6.1) in the class C 2 f _ 2 -2k_ (M E ) ® W(M £ ). We canonically decompose FasF = WQW 1 -, where 
W 1 - represent the orthogonal complement of W in W. Consequentely, u admits the decomposition 

u = u + u T + u , 

where u + u T £ C 2 '^ „_ 2k (M £ ) © W(M e ), whereas u 1 - £ W^. The aim is thus to find a correction 

2h~ 

z e C -x _ ® such tnat > the function to defined by 

°'< 2k 

w := u + u T + z, 
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is a solution to (6.1). It is clear that, as soon as we are able to produce such a correction, then the 



surjectivity in the smaller space C 



2,13 
(5,7- 



(M e ) © W(M e ) will be achieved. Using the linearity of the 



problem and the fact that u is already a solution, one can easily deduce that the function z must satisfy 

L(u e , <?) [z] = h(u E ,g)[u ] in M e . 

A remarkable feature of problem above is that the right hand side is supported by construction only in 
■Di7i,±iJi U D n2y ±R 2 . More precisely, recalling that R[ > Ri and R' 2 > R2, we have that h(u e ,g) [u ] is 
actually supported in the four annuli [R[ - 1, R[) x S n_1 , [--Ri, -iZi + 1] X S™" 1 , [i^ - 1, R' 2 ] x S"- 1 and 
[— i?2, — i?2 + 1] x As a consequence, we are led to solve the following kind of problem 



(6.7) 



L(w e ,g) [w -1 ] in D VuRl , 
on dD VliRl , 



withui,+ e Clf (An^J 8 W+(D ni!R ). Analo gous problems should be posed, respectively, on D m ^ Rl , 
D mt R 2 , D V2t ^R 2 , leading to the construction of vi t -,V2,+,V2,-- Problem (6.7) is, modulo the use of the 
conformal cquivariancc property (2.4), of the same kind of the Dirichlct problem (5.21). Thus, Proposition 
5.6 applies giving the well posedness of (6.7). Once we have obtained these local solutions, we define the 
candidate correction z as 



Vi- + V\- 

vc 

"2,+ + V 2 ,+ 
V2.- + V2,- 



in 
in 
in 
in 
in 



C e := M e \ (D VltRl U D Vu 

Dq 2j — R 2 , 



U D 



U D 



»)2,-fl2 J 



L(u £ ,g) = 

«1,+ = 01,+ 



<9Aji,-Ri 



where i>i,+, i>i,_ , V2.+ , u 2,— 5 w c are the solutions of the following problems 
(6.9) 

(analogous problems for ^2,+, ^2.-) and 



(6.10) 



The Dirichlet data for the problems (6.9) and (6.10) must be chosen in a proper way. Namely, the 
choice of "0 := ( 01,+, tpi,-, 02, +, 02,-) is dictated by the requirement that the function z has the correct 
regularity for being a solution on the whole M e . In fact z is certainly continuous for any choice of 0, but 
there may be a lack C 1 -regularity across the interface dD Vl>Rl . To avoid this situation, we impose the 
following C 1 -matching condition 



7) i v c] 


= 


in 




v c = 


0i,+ 


on 




vc = 


0i,- 


on 




vc = 


■02,+ 


on 


dD m ^ R2 


vc = 


02,- 


on 


dD m _R 2 



.11) 



d ri (vi,± +wi,±)| ri=±fil = d ri v Ciri=±Rl 

dr 2 (v 2 ,± + V 2 ,±)\r 2 = ±R 2 = d r2 V C ^ 2=±R2 - 



To show that the ansatz above actually yields a good definition for z, we adopt the following strategy. 
First of all, we show that problems (6.9) and (6.10) have a unique solution with e-uniform a priori estimate 
for generic Dirichlet data. Secondly, studying the behavior of the so called Dirichlet to Neumann map, 
we will calibrate the choice of the boundary data in such a way that conditions (6.11) are satificd. 
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As we did for (6.7), we note that, modulo the use of the conformal cquivariance property, (6.9) is of the 
same type as (5.37). Thus, Proposition 5.8 applies giving, for any Dirichlet datum ?/>i i+ G C 2 '' 9 (dD VltRl ), 
a unique G C 2 ^ (D Vl;Rl ) © W + (D VuRl ) such that 

( 6 - 12 ) ll S l-+llcif(^ llRl )ffi>V+(^ 1 , Rl ) ^ C W^+\\c^(OD vl , Rl )- 

Of course, the same holds for Vi t -,V2,+, and «2,-- 

To solve problem (6.10) we take advantage of the fact that for any fixed e, the domain C e is compact 
and thus we can apply the Fredholm alternative. Hence, we have existence and uniqueness for problem 
(6.10), provided 

j h(u E ,g)[v]=0 C e 
j v = dC e 

only admits the trivial solution. To prove this, we are going to show that there exists a positive constant 
B > independent of e such that the a priori bound 

(6.13) IMIc 2 -V 2Jt (c e ) - B \\^c° L \ [c e ) 

1 2fc 

is in force for solutions to 

(Uu £ ,g)[v} = f C £ 
\ v = dC e 

To prove (6.13), we use a blow-up argument similar to the one used in the proof of Proposition 6.2. The 
only difference is in the treatment of the case 1. In particular, following the argument used in the above 
mentioned proof and one ends up with a function such that u^ 1 is a non trivial smooth solution 
of the following boundary value problem 

fL(l, 5l )[ii^oo] = in (-Rx, R ± ) x S"- 1 =: C ,i , 
v = on dCo,i 

Thus, thanks to the conformal equivariance property (2.4) and Proposition 5.9, we infer that Voo = 0, 
which is a contradiction. Thus, there exists a unique vc solving (6.10) and such that 

(6.14) || 

Vc 'Wc 2 ' l3 n _ 2k (C E ) — B \\' t P\\c 2 -f 3 (dC E ) ■ 

Having at hand the functions £>i.-, «2.+ , ^i- an d «Ci w e can introduce the Dirichlet to Neu- 
mann maps for problems (6.9) and (6.10). For the exterior problem (6.9), we define the mapping 
T : C 2 -P(dC £ ) — > C x '^{dC e ), whose action is given by 

T : $ := (i/j lt+ ,tpi-, V>2.+ ,-02 ,-) ' — > (-d ri Vi >+ ,d ri v 1 -)-d r2 V2,+,d r2 V2,-)\ +R ■ 

We notice that the action of T is diagonal in the sense that the Dirichlet datum defined on a con- 
nected component of the boundary dC E is mapped to a Neumann datum defined on the same connected 
component. 

On the other hand, for the interior problem (6.10) we define the mapping S e : C 2, P(dC e ) — > C 1 ^(dC e ) 
which acts in the following way 

S e :^i — > (-d ri vc,d ri v c ;-d r2 vc,d r2 vc)\ +R , • 

\r 1 =±R 1 ;r 2 = ±H2 

In terms of the operators T and S e the C 1 -matching condition can be rewritten as 

[(T-S c )[j>]] lt± = ±d Tl v l<± ondD Vu±Rl 
[(T - S £ )[4>}] 2 ± = ±d r2 v 2 ,± ondD n2y±R2 . 
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Thus, the ansatz for z actually yields a well defined correction if we prove that the above pscudodiffcrcntial 
problems is solvable. To this end, we are going to show that, up to choose the parameter e small enough, 
the map 

T-S e : C 2J3 {dC e ) — > C^{dC e ) 

is invertible. First of all, we notice that T and S £ are well defined and, thanks to the a priori estimates 
on the solutions to problems (6.9) and (6.10), they are also e- uniformly bounded. We prove now the 
following 

Lemma 6.6. As e — > the operator S e converges in norm to the operator Sq defined as 
S : C 2 ^(dC ,i) x C 2 ^{dC^ 2 ) — > C^(dC 0A ) x C l 'P(dC , 3 ) 

(("01, + , 01 ,-), (02,+ , 02 ,-)) ' > ((-d ri Vc,l,d ri V C ,l)\ i=±Hi , {-dr 2 VC,2,d r2 V C ,2)\ ±R2 ), 

where the function vc,i is the unique solution of 

!L(l,.gi)[u^ 1 w c .i] =0 in C ,i , 
u Cil = if, li+ on {R x } x g™" 1 
u Cil = ^ _ on {-Ri} x S™" 1 
Of course, an analogous problem characterise vc.2- 

Proof. The proof is by contradiction. Taking advantage of the uniform a priori bound (6.13), we can 
suppose that if the statement is not true, then there exist a sequence (Ej, tpj) such that Sj — ¥ and, for 
every j e N, \\i>j\\c^{dC Bj ) = 1 and 

(6-16) \\(S Sj -S )(^)\\c^(ac e ) > 1/2. 

Let xP c be the unique solution to problem (6.10) with 0j as Dirichlet boundary datum. Up to choose 
a subsequence, we may suppose that the boundary data converge to some in C 2 (dCo y i U dCo^)- 
From the uniform a priori estimates (6.13) combined with the fact that g £j — > gi in C 2 on the compact 
subsets of D m \ {p^}, we deduce that, up to a subsequence, also the functions Vq converge to some 
functions Uqo^ on the compact subsets of Co,i \ {pi} with respect to the C 2 -topology, i = 1,2. Making use 
of the conformal cquivariancc property combined with the removable singularities Lemma 6.1 (which is 
in force since 7 < (n—2k)/k and 2 < 2k < n), it is not difficult to see that Voo,i can be extended through 
Pi to a smooth solution of problem (6.15) on the whole Co,i, i = 1,2. Since it is evident that problem 
(6.15) has a unique solution, the functions Uoo.i an d W00.2 must coincide with vc,i and vc,2 respectively. 
As a consequence, their normal derivative at the boundary must agree. This contradicts (6.16). □ 

In force of the lemma above, the invertibility of T — S e is now a consequence of the invertibility of the 
limit pseudodifferential operator T — Sq. Now, since the spectrum of the limit operator T — Sq is discrete 
it is sufficient to prove the injectivity of (T — So). We reason by contradiction and we assume to have 
i\) = (0i.+ , 02.-, 02.+, 02.-) 7^ for which (T — S*o)[0] = 0. The existence of such a boundary datum 
implies the existence of a smooth function v solving 

L(l, <7i) [w^w] = on D, n , 

and such that 

v(R u 9) = 
v(-R u 0) = 0i,+ (0). 

Moreover, v has the following form 

!Vi i+ in D VuRl 

vc,i in C ,i 

vi- in D Vu - Rl , 
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where we recall that v x,+ and vx,- are solutions to problems of the type (6.9) and belong to C 2 2g (D^^ ) © 
W + (D Vi .r 1 ) and to C^g '(D r , lt -n 1 ) © W~(D 7)1i -r 1 ), respectively. Of course v has the corresponding 
features on D m , but since the situation is somehow symmetric we just focus on D Vl . To reach the desired 
contradiction, we aim to show that v = 0. 

We perform the usual separation of variable, projecting the equation along the eigenfunction <^-'s of the 
Laplacc-Bcltrami operator on (S 1 ™ -1 ,. gg-n-i), having in mind that the high frequencies (j > n + 1) and 
the low frequencies (j = 0, . . . , n) will be treated separately 

n oo 

v(n,e) = ]T^(n)0#)+ V{rx)4>M- 

J=0 j=n+l 

The high frequencies are easier to treat. In fact the deficiency space components are not present in this 
regime. Thus the v 3 are exponentially decreasing for |t"i| — > +oo. Hence, the maximum principle forces 
them to be identically zero. 

To obtain the same result for the low frequencies, we note that the functions v 3 'for j = 1, . . . ,n can be 
written as a linear combination of (ri, 0) and <M'~(ri, 6). Namely, for any j = 1, . . . ,n there exist 
real numbers Aj , Bj such that 

t?'(n) := Aj [*=§kv D , Vl {rx) + VD.Mje^+B, [B=™v D , m (n) - Vd^M)] ^ ■ 

Now, since v 3 = v{ + for r\ > R\ and v\ : + has a component decaying like e~ Sri with S > 1 and the other 
one decaying like e~ ri as n — > +oo, there holds that necessarily Aj = 0. The same argument used when 
r\ < —R\ shows that also Bj = 0. This implies that %p{ ± = for j = 1, . . . , n. 

The last case is when j = 0. As before, the function v° is a linear combination of the two Jacobi fields 
<f>^~ and < &^ + , namely, there exist real numbers A$ and i?o such that 

v^irx) = A <P° n +(n)+ B t>° n -( ri ) . 

Comparing the asymptotic behavior at ±oo of the expression above with the one prescribed by the 
constraints v° — w° + for r\ > R\ and v° — _ for r\ < R\ , we get Aq = = Bq. As a consequence 
it = and V± = 0.' 

The conclusion is that T — So is injective, hence invertible and for e sufficiently small also T — S £ is 
invertible. As already observed, this implies that the correction z is well defined and thus the operator 
h(u e ,g) is surjective on C 2 f _ I ^k(M e ) © W(M e ). 

2k 

To conclude the proof of our statement, we need to establish the a priori estimate (6.6), which obviously 
implies injectivity. First of all we notice that thanks to (5.37), the solution w verifies, for some positive 
constant C > independent of e, the bound 

II^Hcif(^ 1 , Rl )ffiW+(^ 1 ,H 1 ) ^ c [l^llc^_ ( „_ 2fc) (M e ) + \\w\acjc^(dc e )] 

with analogous estimates on the other connected components of M s \ C e . Moreover, the same argument 
used to prove (6.14) implies that for e sufficiently small 

^c 2 -^_ 2k (c e ) ^ c [ Wf\\c°jl _ (n _ aW (iW.) + W W \dC e \\c^(dC e ) ] , 

~< 2fc 

for some C > independent of e. Finally, it follows from standard interior Schauder estimates that the 
norm of w on dC e is uniformly bounded by the norm of / in a small neighborhood of dC £ . Combining 
this remark with the last two estimates, it is now easy to obtain (6.6). This completes the proof. □ 
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7 Nonlinear analysis 

In this last section we are going to correct the approximate solutions u £ to exact solutions, provided the 
parameter e is small enough. Again, for sake of simplicity, we limit oursclf to the case M £ — D ril f l£ D ri2 ^ 
without discussing the minor changes needed for the general case. According to the linear analysis, it 
is natural to look for a correction lying in C 1 ^ n _ 2k (M e ) © W(M £ ). Recalling that W{M e ) is finite 

dimensional and identifying a function in this space with its coordinates with respect to the Jacobi fields 
basis, it is immediate to obtain the following isomorphism of Banach spaces 

J : C^f _„_ 2fc (M e ) © W{M £ ) — > C 2 f _„_ 2fc (M e ) x x R n+1 x M" +1 x M ,l+1 

°lT 2k °>T 2fe 

w = w + a i j ' + Wj+ + a? j '~ i $ij- i — > (w ,a 1 ' + ,a 1 '~,a 2 ' + ,a 2 ~) 
where, for i = 1, 2 and j = 0, . . . , n, 

a j : = XR^j and a j : = X-R'^j and « : = ( a ! ■ • ■ > a n J • 

To describe our perturbation, we denote by u s (a ' + , a 1 ' - , a ,+ , a 2 ' - , •) the variation of u e with param- 
eters a 1 ' , i = l,2, supported on M £ \ C £ , which is defined on D^^, as 

|0 - (a{' + , . . . , 4 ,+ )e" ri ( ri + log |0 - (&}' + , . . . , ^e^ 1 1 + log(l + aj'+) ) . 

The definition of M £ (a 1,+ , a 1 ' - , a 2 ' + , a 2 ' - , • ) on the other connected component of M £ \ C £ is analogous. 
We note en passant that the 'straight' approximate solution u £ (-) corresponds to a l,± = 0, for i = 1,2. 

To get exact solutions for our nonlinear problem, we are led to look for (w, a 1,+ , a 1 '", a 2 ' + , a 2 '~) g 
C^f _„^(Af 6 ) x R n+1 x W l+1 x R n+1 x R" +1 such that 

°lT 2fe 

(7.1) ^(^(a 1 ^, a 1 '" o 2 '+a a '- •)+«>(•)» fl) = °- 

A simple computation gives 

A/" (w e (sd 1,+ , set ,_ , sd 2,+ , sd 2: ~, ■) + sw(-),g) — L(u e (-), g)[w] , 



d 
ds 



s=0 

where, according to (7.1), 

(7.2) w = w + a^ + ^;+ + a, tf* . 

This formula suggests that the perturbation of u £ to an exact solution will involve a decaying term 
(through w) together with small conformal variations of the former 'straight' approximate solution. Geo- 
metrically, these variations corresponds to translations along the Delaunay axis, changes of the Delaunay 
parameter rj and 'bendings at infinity'. 

Using a Taylor expansion, we can rewrite (7.1) as 

= Af(u s (a 1>+ ,a 1 '-,a 2>+ ,a 2 '-,-)+w(-),g) 

= M(u £ (-),g) + h(u e (-),g)[w] + Q(u e (-),g) [w,w] , 

where Q(u e (-), g) [w, w] is the quadratic remainder. Thus, the fully nonlinear problem becomes equivalent 
to the following fixed point problem for w = (w, a^+, a 1 ~, a 2 >+, a 2 ~) g C 2 f (M £ ) xW n+1 xl" +1 x 

(7.3) w = L(u e (-),ff)- 1 [ -M(u £ (-),g) - Q(u £ (-),g)[w,w 



2G 



It is worth remarking that at first time one could have used a simple perturbation of the form u £ + w 
with w as (7.2). In fact, the first order expansion of N(u £ + w,g) is also given by L(i/ £ (-), g)[w] and 
thus the linear analysis of the previous sections is still in force. On the other hand the components of 
the correction w along the Jacobi fields are not necessarily infinitesimal with respect to u e when 

\ri\ — > +00, i = 1,2, and this may possibly affect the completeness of the final solution. To avoid this 
problem we had to deal with perturbations of the form (7.1). Indeed the conformal equivariance of the 
exfc-equation insures that the variations u e (a 1 ' + , a 1,_ , a 2,+ , a 2 ' - , • ) are still complete exact solutions far 
away from the central body C e . Since the remaining part of the perturbation w is exponentially decaying 
at infinity, the completeness of the exact solutions is definitely guaranteed. 

We denote by V the mapping V : C 2 f _„_ 2fc (M e ) x R" +1 x R n+1 x M ,l+1 x M" +1 -> _ n _ 2k (M e ) x 

R n+1 x R n+1 x R' i+1 x W l+1 that associates to any to the right hand side of (7.3). In what follows, we 
will find the fixed point w as a limit of the sequence {tfi}i6N defined by means of the following Newton 
iteration scheme 



(7.4) 






V(wi), i G N. 



To this end, we need some preparatory Lemmata. We state the following 

Lemma 7.1. There exists a positive constant A = A(n,k) > such that for every 1 < S < 8{n,k) 
and < 7 < (n — 2k) /k the proper error is estimated as 



(7.5) \\M(u £ ,g)\\ c o., , < Ae^ 2 ' 

— 0.7 — (m — 2fc) ^ fc ' 



The proof of this result could be found in [4] to which we refer for all the details. Incidentally, we notice 
that the proof substantially uses that Af(u e ,g) is concentrated only on the neck region N £ . Thus, even if 
we have to deal with noncompact manifolds, the computations needed to estimate this term are exactly 
of the same type of the computations in [4, Lemma 5.1]. 

In the following lemma we provide an estimate for the quadrate remainder outside the neck region. 

Lemma 7.2. There exists a positive constant C = C(n, k, 6, 7) such that for every 1 < 8 < S (n, k) 
and < 7 < (n — 2k) /k there holds 

(7.6) \\w\\ C 2,f> _ 2fe(Me)eW(ME) <P =>■ \\Q(u E ,g)[w,w}\\ c oj^_^_^ {MANs) < Cp 2 . 

5,7 2k 

Proof. Let us fix a positive p for which 

IMIc 2 -^ _ s .(M £ )eW(M 6 ) ^ P 

- S -1- !!L 2k— 

holds. We analyze the size of Q(u e ,g)[w; w \ according to the definition of the norm in C^'J 1 , n _ 2fc j (Af e ) 
in (6.2). In particular, we decompose M e \N e as M e \N e = Ufl I|lr H I UD^r^ UD W) _jj, UC e \N E 

and we prove that (7.6) holds on every component of the above decomposition. On this regard, let us 
notice that it will be sufficient to check (7.6) only on C e \ (N £ n D Vl ) and on D m n 1 . 
We start with the analysis on D^r^. Wc recall that in this region g — v ^/^ n ~ 2k ^ g cyl _ Hence, from the 
computational point of view it is more convenient to work with the cylindrical metric as a background 
metric. To this end we set z := vjj^ w, z := vjj^ w and z T := »n w T . Thus, by using the conformal 
equivariance property (2.4) and that u s = 1 on D m _R 1 (see (4.1)), one has 

(7.7) Q(u £ (-),g) [w;w] = v D ^ 2k Q{v Dlii {-),g cy i) [z; z], w = v^z. 
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Since vd v1 is uniformly bounded from above and from below, one can deduce the desired estimate from 
its analogous for Q(vn m {),gcyi) [z ; z\. This last term can be expanded on D m ji 1 as 



Q(v Dm (■), 9c V i) [z;z] = Q(v Dm (•), g cyl ) [(z, a 1 ^, 0, 0, 0) ; (z, 0, 0, 0)] 

h(v Dm (a 1 ' 4 ", 0, 0, 0, •) + rz(-),g cyl ) - h(v Dtll (a 1 ' 4 ", 0, 0, 0, -),g cy i) 

H v D m (a 1, + , 0, 0, 0, •), 9c V l) - Mv Dvi (•), gcyl) 

(7-8) +/ [L(«zj, 1 (TO 1 ' + ,0 I 0,0 J -) J ffc»i)-M«o, 1 (-).ffc»i) [a]' + 9^]dr 



(It 



To proceed, we recall that the linearization of Af(-,g C yl) around vo ril has the following general structure 
(see (2.1) and [17, Section 5]) 



(7.9) 



H V D m , gcyl ) = L° (VZV > gcyl) + C^kVB^" 



where c n ^ is a computable positive constant and L°(v£> j;i , g cy i) is a second order differential operator 
with smooth coefficients of the following form 



(7.10) 



^°(v Dvi ,g cyl ) = P^ivD^Vvn^VnJd", 

\a\<2 



where a is a multi-index and P 2k 1 : I 
2k- 1 

Pl k -\x,y,z) :-- 



i8o+|/9 1 |+| J 8 3 |=2fc-l 



As a consequence, setting h := (ho, hi, h%) elx M n x 
has 



is an homogeneous polynomial of degree 



and expanding at first order P 2 one 



(7.11) Pa(x + h ,y + h 1 , z + h 2 )-P^- L (x,y,z) = DPf- L (x,y,z)-h + 0(\h\'). 

We have now all the ingredients to obtain the estimate on D m ^ 1 for (7.8). First of all, we recall that 
(see (5.18) and (5.19)) 

\\Q{vD vl ,9cyi)[z;z]\\oM {D ) := sup (coshn) 5 \Q(v D ,g cy i)[z;z}\ 



(7.12) 



+ SUP (cOShri) [Q(v D ,g c yl)[z]z]] C 0,f},, r ^_ lri+1)xB n-ly 

ri>ili+l 



We will estimate separately the two terms in (7.12). 

We start with the estimate of the weighted C° norm of Q\. To this end, by applying (7.11) to the 
operator L° and simply expanding at first order the remaining term in (7.9) we may decompose Q\ into 
Qi = +91,2 where, for any (r±,6) G D vltRl , 

(7-13) <z M := f J2 [DPl k -\v D ^Vv D ^ 2 v D J-Th + 0{\Th\ 2 )}d a zdT 



l«|<2 



and 
(7.14) 



91,2 := / (d n . k v^ k tz + 0(\tz\ 2 )[z] dr, d n , k := c„, fc - l) 



28 



where the vector h appearing in (7.13) has components (z, Vz, V 2 z). Thus, it is immediate to obtain 

(7.15) ||gi,l|| C ° 5 (A/ e \C e ) < C \\z\\ 2 Cl s (M c \C c ) and hia\\c°_ s {M B \Ce) < C Pfc2_ s (M e \C e ) 

for some positive constant C, possibly depending on n, k, 7 and <5. 

Concerning Q2, we preliminarly expand (a 1,+ , 0, 0, 0, • ) as (recall (5.15)-(5.17) and (5.35)) 

(7.16) v Dvi {a}' + , 0,0,0,-) = w^JO+^o^ + Odl^ll^.)). 
Thus, splitting Q 2 into Q 2 = 92,1 + 92,2, where 

(7.17) fe! := [L°(^ i (-) + ^+a]- + + 0(||z T ||2 v( ^ ) ), 5c , ( ) - L°(^ {-),g cyl ) ] [z] 
and 

(7.18) 92,2 := C n, fc [(^J0 + n + ^ + +°(^ T |lW)))^^- c ^ v Sf" _1 (-)] W 
Thus, by using (7.11) in (7.17) and expanding at first order in (7.18) we get 

||92,l||c° 5 {M e \C s ) < C\\z T \\ w{Me) \\z\\ c 2jj ( Me )ffiW(M e ) 

-"'1 27T" 

(7-19) \\<l2,2\\c° s (M E \C e ) < ^IMIw(M,)|W| C ^ n _ 2k{Mc)(SW{Me) , 

where the positive constant C possibly depends on n, k, 7 and S. 

Now, we estimate Q3. The estimate relies on the observation that Q3 has compact support. To see this, 
let us show that 

Q'3 ■■= [U«i?,Jo 1,+ .0,0 > O r ), flctf O-M«iJ, 1 (0 ) ffcyO][5 J - ,+ *i + ] 
has indeed compact support on D m ji 1 . We may decompose Q' 3 as 

Q'a = [M«2J, 1 (o 1 ' + > 0,0 > 0,0,5cyO-M«D, 1 (o 1,+ ) > 0,0 > -) > ff OI ,{)] [S) ,+ *^ + ] 
+ [ L(« Dni (a^+, 0, 0, 0, •), 9c V l) ~ Uvd v1 (•), 9eyl) ] [(a/ + - <# + )tt*+] 
+ [Uv Diii (a 1 '+,0 7 0,0r),9cyi)-Uv D , ii (-) 7 g cyl )] [a}' + **+] 

= : 93,1 + 93,2 + 93,3 • 

Now, q 3 3 = in D rilt n 1 . In fact, we observe that 

9 3 ,3 = / £ 2 ^(«D„>a 1,+ , 0,0, 0,-),^) [a^4+,a l 1 ' + ^+] 





and that for any (n, 0) S D m ^ 1 

M{v Drii {a}> + ,Q,Q,Q,ri,6),g cyl ) = 0. 

Now, since aj' + — a*' + = aJ'^Xi^ — 1) has compact support, it turns out that also q' 32 is compactly 
supported. To see that the remaining term q' 3 1 has compact support, we first expand VD m (a 1,+ , 0, 0, 0, •) 
around Vd v1 (a 1,+ , 0, 0, 0, ■) as 

v Dni (^+, 0, 0, 0, •) = v Dm (a 1 -*, 0, 0, 0, •) + - a[ + ) + Q(\a^ a 1 ^ 2 ). 
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As before, note that dv ° vi ^° a i,+°'°' 0, - ~ a ]' + ) + 0(|a 1,+ — a 1,+ | 2 ) has compact support. Now, 

recalling (7.10) and (7.11) and expanding at first order the potential term in (7.9), it is not difficult to 
get 

3 

which clearly implies that q' 3 1 is compactly suported. 

We can now give the desired estimate for Q3. In particular, thanks to the above computations, it is 
evident that we can equivalcntly estimate the C° norm of Q3 instead of its C°_ s norm. To obtain this 
estimate we reason as before. Using (7.16), (7.10) and (7.11) and expanding at first order the potential 
term in (7.9) it is standard to get 

(7.20) sup \Q 3 \(n,e) < C\\z t \\ 2 w{Me) . 

(n,0) e (i?.i,+oo)xS™-! 

Thus, collecting (7.15), (7.19), (7.20) and recalling that Dp is uniformly bounded from below and from 
above, we get (see (7.7)) the weighted C° estimate for Q(u e ,g)[w;w] on D Vi r 1 , namely the following 

(7.21) sup (coshn) 6 \Q(u e ,g)[w;w]\(n, 6) < Cp 2 . 

Now, we turn our attention to the estimate for the Holder quotients. We will use two different strategies. 
In particular, for the terms Q\ and Q2, we will estimate directly their Holder quotient. For Q3 we will 
estimate its weighted C 1 norm. This is possible thanks to its particular structure. More precisely, it is 
possible to obtain a weighted C 1 estimate by relying, loosely speaking, on the regularity of the Jacobi 
fields and of the a}' + := XR' al,+ - 

We start with the estimate of the term Q3. By first using (7.16) and then expanding at first order the 
coefficients of the linearized operator as in (7.11), it is sufficient to get a weighted C° estimate for 

V [ J2 [DPl k -\v Dn ,Vv Dn ,V 2 v Dn )-h+0{\h\ 2 )} d a {a\ + ^+)} 

M<2 

and for 

V [ {d n , k v D ^(-)^~ 2 h +O(\h \ 2 ) [a]' + ^] ] , d n , k := c„ ifc (^ - l), 

where the vector h has components hi = V l (a J -' + ^' 1 + + 0(||w T ||>v(.M ))' ^ or * = 1, 2. We will outline 
only the estimate for the first term. A similar argument applies to the second. First of all, from the 
definition of d 1,+ := XR' l al ' + (recall that the cut off function XRi is smooth and bounded with its 
derivatives) and from the definition of the Jacobi fields, we easily get 

|V^(a}'+^>+)| < C||z T || w(A4e) , \a\ <2 
Thus, since sup ri>fll |V 1 u_d 7ji | < C, for i = 0, 1, 2, 3, we get 

sup |VQ 3 |(ri,0) < C\\z t \\ 2 w{Me) , 

which implies, together with (7.20), 

(7-22) WQ»\\c°ji(D n ) < Cll^llwt^)- 

Now, we estimate the Holder quotients for Q\ and Q^- We start with Q\. As we did for the C° estimate, 
we split Qi into Qi = 5x1 + q% 2- We will detail only the estimate for the Holder quotient of g_i\, the 
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dr, 



one for being completely analogous. Recalling that h is the vector with components hi = Vg c z, for 
i = 0, 1, 2, we can write 

70 H<2 

- DPf (r'), V« DfI1 (r'),V 2 v Dvi (/)) [rh(r', 0')] d Q i(r', 9') 

{7.23) + f V [G(r/i(r,6»))a Q z(r,6l)-G(T/i(r',6l , ))<9 Q z(r',6l') 

Jo |a|<2 L 

where G is a smooth function such that G(v) = 0(\v\ 2 ) and DG(v) = 0(\v\). Using the short notation 

A a (r,6) := DP 2k -\v Dlll (r), (r), V 2 ^ (r)) , 

we split integrand of the first summand in the expression above into 

£ [A Q (r,0)[r/i(r,0)]-[9 Q z(r,0) -3 Q z(r',0')] + A Q (r', 9') [rh{r, 9) - rh(r', 9')] ■ d a z{r', 9') 

+ [A a (r, 9) - A a (r', 9')] [rh{r, 9)] ■ d a z(r',9') 

Using the fact that (coshr)~ A < 1, for 8 > it is now easy to bound the weighted Holder quotient of 
each term by a constant times llzll 2 2 » . Applying the same reasoning to the second summand in 

C -S ( D vi-Ri) 

(7.23) and to the term qi^, one concludes that 

(7-24) sup (coshr 1 ) i [Q 1 ] c0 ,, (( 1)xS „_ 1} < Cp|j^ 2 ,, 

r 1 >B. 1 + l U -J 

Using the same arguments, one can deduce the same type of estimate for the Holder quotient of Q2, 
namely 

< ~ 7 ' 25 ' ) r XR P +1 ( C0Shri ^ [ < 5 2 '] C0 '^('- 1 ' r + 1 ) ><S " _1 ) - C \\2\\c™ (M E \cJ\ Z Wc 2 ^(M E \C e )(BW(M E )- 

Thus, combining (7.21) with (7.22), (7.24) and (7.25) and recalling that ud is uniformly bounded from 
above and from below, we obtain 

(7.26) \\Q(u e ,g)[w;w]\\ c o,p iDt)iRi) < Cp 2 . 

As anticipated, the estimates of Q(u E (-),g)[- ; •] on the other ends D m ^-n 1 , D mj R 2 , D n2t -R 2 clearly 
follows from a similar argument. Thus, (7.26) actually becomes 

(7.27) IIQ(^,s)K™]llc°f(M s \c e ) < °P 2 - 

Finally, it remains to estimate Q(u E (-) : g) on C E \ N E . Since on this region u E (a 1 ' + , a '~, a 2 ' + , a 2 ' - , ■ ) 
coincides with u E (-), it turns out that the quadratic remainder can be written as 

Q( u e(-),g) [w, w] = Q(u E (-),g)[w;w) 

r-l 

[h{u E {-) + tw(-), g) - h(u E (-),g) ] [w] dr . 

Thus, using an argument similar to the one used above (alternatively, one may refer to [4]), we have 

\\Q(u E ,g)[w;w]\\ c o,e {Ce \ Ne) < Cp 2 . 
Thus, the lemma is proven. □ 
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We are now in the position to conclude the proof of Theorem 1. We need to prove that the sequence 
of the solutions to the iterative scheme (7.4) (which exist thanks to Proposition 6.6) is equibounded in 
C -x _^2h( M e) © W(M e ). We start with the estimate on w\ . Thanks to the uniform a priori estimate 

(6.6) for the linearized problem and to the estimate of the proper error term in Lemma 7.1, we immediately 
have 

(7.28) l|wi|| c 2 ^ „_ 2fe (M e )ew(M e ) < ALe^ 2 ^, 

where the constant L = L(<5, 7,71, k) denotes the uniform bound on the norm of L(m £ (-), g) , while the 
constant A = A(5, 7, n, k) is the constant appearing in Lemma 7.1. 

We proceed with the estimate of W2- From the very definition of W2, we have 

(7.29) \\w 2 \\ C 2,f> (Me)eW(Me) < L\\M(u e ,g) + Qiue^^w^wi)]]^ {Me) 

S '~< 2k ' 

< ALe^ +2 ^ + L||e(u 6 ,fl)[«; 1 ;«;i]|| c o, fl , M y 

— 5,7 — (n — 2fc) ^ e ' 

Thus, we need to estimate the quadratic remainder. Recalling the definition of the global weighted norm 
(6.2) in M e , we have the following 

\\Q(u e ,g)[wi;wi]\\ c o,p (M s = \\Q(u e ,g)[wi;wi]\\ c o,p (M XN ) + sup N (e coshi) 7- ^ 1-2 ' ) I Q(u £ , g)[wi; wi] 

+ su Pte(io ge ,-io ge ) (ecosht)T-("- 2fc ) [ Q{u e ,g)[w x ] ] c w{ ( t -i,t+i)xs^ ) 
Thanks to Lemma 7.2 and to estimate (7.28), the first term is estimated in this way 

|| Q(u E ,g)[w 1 ;w 1 }\\ c o, !{MeXNe) < C\\ Wl \\l^ {Me)mW(Ms) < ACLs^+ 2 ^. 

We consider now the second term (the term containing the Holder quotients will be then estimated in 
the same way, see [4]). On the neck region N £ , Q(u £ (-),g)[wi;wi] has this form 

Q(ue(-),S) [wi;wi] = Q(u e (-),g)[u> 1 - ) u) 1 ] 

[L(Ue(.) + TW(-),g) -L(Ue(.))][t/)l]dT. 
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Thus, by exploiting the structure of L(it e (-), g) given in (7.9), using a first order expansion and recalling 
also (7.28) we have that there exists a positive constant independent of e such that 

(£cosh^- ( "- 2fc > \Q(u E (-),g)[wx;wi\\ < C(ecosht) 7 - ( "- 2 ^(ecosht) (2fc ^ 2 ) I1 ^(£coshO _27+ ^l|w'i|| 



(7.30) < ACLe-^e^ +2) - 



-6,-y- 



,. (M e )eW(M e )' 



where we have used also the fact that, for j = 0,1,2, V^(it e ) = 0((e coshi) "21- ) on the neck region. 
Now, since —7 + (7 + 2) "~ 2fc > 0, for any 7 e (0, n ~ 2k ), we have that, setting 

B := AL 2 Ce-^s^+ 2 ^, 

there exists a positive number £0 = £o(S,^,n,k) such that, for any e £ (0,£o], we can choose B < \. 
Consequentely, the estimate (7.29) for w 2 becomes 

I|W2|| C 2 ^ n _ 2k (U e )mW(M e ) ^ L\\M(u E ,g) + Q(u E ,g)[ Wl ; Wl }\\ c o^^^ {Me) 



< ALe^+ 2 ^ + L\\ e(« 8 ,5)K;«;i]|| o,/i 



a , T - (re -2*,)( M =) 
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Then, we can iterate the above estimate obtaining, for j > 1, 



z~2k 



(7.31) \\w j+l \\ G 2,,^ n _^ {Me)BW{Me) < ALe^+V 1 » o i+lj 

where the sequence dj is inductively defined as 




1 + Ul JGN. 



1 

Now, a straightforward induction argument shows that sup^ cij < 2, thus estimate (7.31) becomes 
(7-32) ll^llc^_ n _ 2fe (M e )eW(M e ) < 2AL^ +2 >^ . 

The previous estimate, combined with the fact that the embedding C 2 '^ (M e ) — > C^_ S ,(M S ) is compact 
for any 5' < 8 (see [18, Chapter 12]), implies (up to a subsequence) the convergence in C 2 ^, (M e ) of Wi 
to a fixed point w e = (w £ , a £,1 - + , a £,1 '~, a e,2,+ , a £ ' 2 '~) of the problem (7.3). 

Thanks to the canonical identification (7.1) wc will write, with a little abuse of notation, w £ — w £ + 
o^' t '+\^j,+ 4- a^ ,t,_ ^'-. Since (7.32) is uniform with respect to j, w £ verifies 

(7.33) |K||c^ tt (M t )ew(M t ) < 2Aie^+ 2 )^ . 

We claim now that there exists eq > 0, such that for every e € [0, eo) the exact solutions 

u e (a e>1 ' + ,o e ' 1 -,o e,2 ' + ,o E ' 2 -,-) + w £ (-) 

are positive. To see this fact we first observe that up to choose e sufficiently small, the function y £ := 
M e (a e,1 ' + , a £,1, ~, a e ' 2,+ , a £,2, ~, ■) is positive everywhere by definition. Secondly, since w £ decays faster 
than y £ along the complete ends, the exact solution y £ + w £ is certainly positive outside of a compact 
region Kq C M £ . Hence, since (7.33) implies that 

(7-34) Wvl'weWc^Ko) < Ce-^+V 2 ^, 

there holds that u £ {a £ ' l - + , a £,1 '~, a £,2 ' + , a £ ' 2 '~, •) + w £ = y £ + w £ = y e (l + y~ x w £ ) > 0. 

For e G (0, eo] we set 

g £ := (u e (a e > 1 '+,a s ' 1 >-,a E > 2 > + ,a E > 2 >-,-)+w e )^g. 

The above considerations imply that g £ is the metric sought. We recall that the completeness of these met- 
rics is a consequence of the decaying of w £ on the ends of M £ and of the fact that u £ (a £,1 ' + , a e,1 '~ , a e,2 ' + , a £ ' 2 '~ , 
is, by construction, a complete solution to the cr^-cquation, locally on the end. 

Moreover, the family of metrics g £ converges to the initial metric gi with respect to the C 2 topology 
on every compact subset of D m \ {pi}, for i = 1,2. This is evident on the four ends D r]it ±R i . In 
fact on these regions we have g = g% and (7.33) implies that on every compact subset of D rii ^±p i 
u £ (a £ ' 1 ' + , a 5,1 ' - , a £ > 2 ' + , a £ ' 2 -~, •) + w £ — > u E (-) = 1 in C 2 as e — > 0. To see that g £ — > gi on D Vi \ 
{pi} n C £ , we recall, as before, that u 6 (a e ' 1,+ , a e,1 '~, a £ - 2:+ , a e ' 2, ~, •) = u £ on C £ . Thus, the metric g £ 
could be written as 

9e = (1 + U^We)^ g e . 

Now, since by construction on every compact subset of D Vi \ {pi} n C £ the metric g £ converges to the 
initial metric gi in C 2 as e — > 0, (7.34) implies that also the exact solutions g £ tend to the initial metric 
gi with respect to the C 2 -topology on the compact subsets of D rji \ {pi} n C £ , for i = 1,2, as e — >• 0. 
This concludes the proof of Theorem 1. 
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